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Abstract. This paper studies the construction of geometric integrators for 
nonholonomic systems. We derive the nonholonomic discrete Euler-Lagrange 
equations in a setting which permits to deduce geometric integrators for con- 
tinuous nonholonomic systems (reduced or not). The formalism is given in 
terms of Lie groupoids, specifying a discrete Lagrangian and a constraint sub- 
manifold on it. Additionally, it is necessary to fix a vector subbundle of the Lie 
^ algebroid associated to the Lie groupoid. We also discuss the existence of non- 

'~H holonomic evolution operators in terms of the discrete nonholonomic Legendre 

transformations and in terms of adequate decompositions of the prolongation 
of the Lie groupoid. The characterization of the reversibility of the evolution 
operator and the discrete nonholonomic momentum equation are also consid- 
" ered. Finally, we illustrate with several classical examples the wide range of 

application of the theory (the discrete nonholonomic constrained particle, the 

OSuslov system, the Chaplygin sleigh, the Veselova system, the rolling ball on 
a rotating table and the two wheeled planar mobile robot). 
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1. Introduction 

In the paper of Moser and Veselov [30] dedicated to the complete integrability 
of certain dynamical systems, the authors proposed a discretization of the tangent 
bundle TQ of a configuration space Q replacing it by the product Q x Q, approx- 
imating a tangent vector on Q by a pair of 'close' points (qo;9i)- In this sense, 
the continuous Lagrangian function L : TQ — > R is replaced by a discretization 
: Q X Q — > M. Then, applying a suitable variational principle, it is possible to 
derive the discrete equations of motion. In the regular case, one obtains an evolu- 
tion operator, a map which assigns to each pair {qk-i,qk) a pair (g^, qk+i), sharing 
many properties with the continuous system, in particular, symplecticity, momen- 
tum conservation and a good energy behavior. We refer to [32" for an excellent 
review in discrete Mechanics (on Q x Q) and its numerical implementation. 

On the other hand, in |40l |44] , the authors also considered discrete Lagrangians 
defined on a Lie group G where the evolution operator is given by a diffeomorphism 
of G. 

All the above examples led to A. Weinstein to study discrete mechanics on 
Lie groupoids. A Lie groupoid is a geometric structure that includes as particular 
examples the case of cartesian products Q x Q as well as Lie groups and other 
examples as Atiyah or action Lie groupoids [1^ . In a recent paper [27] , we studied 
discrete Lagrangian and Hamiltonian Mechanics on Lie groupoids, deriving from 
a variational principle the discrete Euler-Lagrange equations. We also introduced 
a symplectic 2-section (which is preserved by the Lagrange evolution operator) 
and defined the Hamiltonian evolution operator, in terms of the discrete Legendre 
transformations, which is a symplectic map with respect to the canonical symplectic 
2-section on the prolongation of the dual of the Lie algebroid of the given groupoid. 
These techniques include as particular cases the classical discrete Euler-Lagrange 
equations, the discrete Euler-Poincare equations (see [H |6l [29j [30] ) and the discrete 
Lagrange-Poincare equations. In fact, the results in |27j may be applied in the 
construction of geometric integrators for continuous Lagrangian systems which are 
invariant under the action of a symmetry Lie group (see also [IH] for the particular 
case when the symmetry Lie group is abelian). 

From the perspective of geometric integration, there are a great interest in intro- 
ducing new geometric techniques for developing numerical integrators since stan- 
dard methods often introduce some spurious effects like dissipation in conservative 
systems [HI [42]. The case of dynamical systems subjected to constraints is also 
of considerable interest. In particular, the case of holonomic constraints is well 
established in the literature of geometric integration, for instance, in simulation of 
molecular dynamics where the constraints may be molecular bond lengths or angles 
and also in multibody dynamics (see [THj [20] and references therein) . 
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By contrast, the construction of geometric integrators for the case of nonholo- 
nomic constraints is less well understood. This type of constraints appears, for 
instance, in mechanical models of convex rigid bodies rolling without sliding on a 
surface [41j . The study of systems with nonholonomic constraints goes back to the 
XIX century. The equations of motion were obtained applying either D'Alembert's 
principle of virtual work or Gauss principle of least constraint. Recently, many 
authors have shown a new interest in that theory and also in its relation to the 
new developments in control theory and robotics using geometric techniques (see, 
for instance, [1 d H |H1 [H HI [H] ) . 

Geometrically, nonholonomic constraints are globally described by a submanifold 
M of the velocity phase space TQ. If M is a vector subbundle of TQ, we are dealing 
with the case of linear constraints and, in the case M is an affine subbundle, we are in 
the case of affine constraints. Lagrange-D'Alembert's or Chetaev's principles allow 
us to determine the set of possible values of the constraint forces only from the set 
of admissible kinematic states, that is, from the constraint manifold M determined 
by the vanishing of the nonholonomic constraints Therefore, assuming that the 
dynamical properties of the system are mathematically described by a Lagrangian 
function L : TQ — > K and by a constraint submanifold M, the equations of motion, 
following Chetaev's principle, are 



Sq' = 



where Sq^ denotes the virtual displacements verifying ——Sq^ = 0. By using the 

oq^ 

Lagrange multiplier rule, we obtain that 

fdL\ _dL - 
dt\dct) dq^~ ""dq^' ^ ' 

with the condition q(t) G M, Xa being the Lagrange multipliers to be determined. 
Recently, J. Cortes et al [9] (see also [HI [38l [39]) proposed a unified framework for 
nonholonomic systems in the Lie algebroid setting that we will use along this paper 
generalizing some previous work for free Lagrangian mechanics on Lie algebroids 
(see, for instance, ^ 1321 |31I35] ) . 

The construction of geometric integrators for Equations ( |1.1[ ) is very recent. In 
fact, in |37j appears as an open problem: 

...The problem for the more general class of non-holonomic con- 
straints is still open, as is the question of the correct analogue 
of symplectic integration for non-holonomically constrained La- 
grangian systems... 

Numerical integrators derived from discrete variational principles have proved their 
adaptability to many situations: collisions, classical field theory, external forces... ^51 
132] and it also seems very adequate for nonholonomic systems, since nonholonomic 
equations of motion come from Holder's variational principle which is not a stan- 
dard variational principle [IJ, but admits an adequate discretization. This is the 
procedure introduced by J. Cortes and S. Martinez [8j [10] and followed by other 
authors [T2J |T31 [TSJ [33] extending, moreover, the results to nonholonomic systems 
defined on Lie groups (see also |25l for a different approach using generating func- 
tions). 

In this paper, we tackle the problem from the unifying point of view of Lie 
groupoids (see [3] for the continuous case). This technique permits to recover all 
the previous methods in the literature [TUl [3 [3S] and consider new cases of great 



4 D. IGLESIAS, J. C. MARRERO, D. MARTfN DE DIEGO, AND E. MARTfNEZ 



importance in nonholonomic dynamics. For instance, using action Lie groupoids, 
we may discretize LR-nonholonomic systems such as the Veselova system or us- 
ing Atiyah Lie groupoids we find discrete versions for the reduced equations of 
nonholonomic systems with symmetry. 

The paper is structured as follows. In section 2 we review some basic results on 
Lie algebroids and Lie groupoids. In particular, we describe the prolongation of a 
Lie groupoid P5| , which has a double structure of Lie groupoid and Lie algebroid. 
Then, we briefly expose the geometric structure of discrete unconstrained mechanics 
on Lie groupoids: Poincare-Cartan sections, Legendre transformations... The main 
results of the paper appear in section 3, where the geometric structure of discrete 
nonholonomic systems on Lie groupoids is considered. In particular, given a discrete 
Lagrangian : F — > M on a Lie groupoid F, a constraint distribution Dc in the 
Lie algebroid Ey of F and a discrete constraint submanifold Mc in F, we obtain 
the nonholonomic discrete Euler-Lagrangc equations from a discrete Generalized 
Holder's principle (see section 3.1). In addition, we characterize the regularity of the 
nonholonomic system in terms of the nonholonomic Legendre transformations and 
decompositions of the prolongation of the Lie groupoid. In the case when the system 
is regular, we can define the nonholonomic evolution operator. An interesting 
situation, studied in in Section 3.4, is that of reversible discrete nonholonomic 
Lagrangian systems, where the Lagrangian and the discrete constraint submanifold 
are invariants with respect to the inversion of the Lie groupoid. The particular 
example of reversible systems in the pair groupoid Q x Q was first studied in . 
We also define the discrete nonholonomic momentum map. In order to give an 
idea of the breadth and flexibility of the proposed formalism, several examples are 
discussed, including their regularity and their reversibility: 

- Discrete holonomic Lagrangian systems on a Lie groupoid, which are a 
generalization of the Shake algorithm for holonomic systems p!51 EUl 152] : 

- Discrete nonholonomic systems on the pair groupoid, recovering the equa- 
tions first considered in [TU]. An explicit example of this situation is the 
discrete nonholonomic constrained particle. 

- Discrete nonholonomic systems on Lie groups, where the equations that 
are obtained are the so-called discrete Euler-Poincare-Suslov equations (see 
[H]). We remark that, although our equations coincide with those in [T¥j . 
the technique developed in this paper is different to the one in that paper. 
Two explicit examples which we describe here are the Suslov system and 
the Chaplygin sleigh. 

- Discrete nonholonomic Lagrangian systems on an action Lie groupoid. 
This example is quite interesting since it allows us to discretize a well- 
known nonholonomic LR-system: the Veselova system (see [H]; see also 
[T^l. For this example, we obtain a discrete system that is not reversible 
and we show that the system is regular in a neighborhood around the 
manifold of units. 

- Discrete nonholonomic Lagrangian systems on an Atiyah Lie groupoid. 
With this example, we are able to discretize reduced systems, in particular, 
we concentrate on the example of the discretization of the equations of 
motion of a rolling ball without sliding on a rotating table with constant 
angular velocity. 

- Discrete Chaplygin systems, which are regular systems (LdjMc'I'c) on 
the Lie groupoid F =J M, for which (a,/?) o ij^^^ : Mc M x M is 
a diffeomorphism and p o : Dc TM is an isomorphism of vector 
bundles, (a, f3) being the source and target of the Lie groupoid F and p 
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being the anchor map of the Lie algebroid E-p- This example includes a 
discretization of the two wheeled planar mobile robot. 

We conclude our paper with future lines of work. 

2. Discrete Unconstrained Lagrangian Systems on Lie Groupoids 

2.1. Lie algebroids. A Lie algebroid E over a manifold M is a real vector bundle 
T : E ^ M together with a Lie bracket |-, •] on the space Sec(r) of the global cross 
sections oi t : E ^ M and a bundle map p : E TM, called the anchor 
map, such that if we also denote by p : Sec(r) X(M) the homomorphism of 
(M)-modules induced by the anchor map then 

lXJYj = flX,Y]+piX){f)Y, (2.1) 

for X,Y e Sec(T) and / e C°°{M) (see [55]). 

If (E'j |-, •], p) is a Lie algebroid over M then the anchor map p : Sec(T) — 
X(M) is a homomorphism between the Lie algebras (Sec(T), |-, •]) and (X(Af ), [•,•]). 
Moreover, one may define the differential rf of _E as follows: 

k 

d^l{Xo, . . . , Xfc) = ^(-i)V(^.)(/^(^o, . . . , ^, . . . , Xk)) 

(2.2) 

+ J2i-ly+'KlX^.XJlXo,...X,■■■,X3,■■■,Xk), 

i<j 

for ^ G Sec(A'^r*) and Xg, . . . ,Xk & Sec(r). d is a cohomology operator, that is, 
~ 0. In particular, \i f : M — > M is a real smooth function then df{X) = p{X)f, 
for X e Sec(T). 

Trivial examples of Lie algebroids are a real Lie algebra of finite dimension (in 
this case, the base space is a single point) and the tangent bundle of a manifold M. 

On the other hand, let {E, |-, •], p) be a Lie algebroid of rank n over a manifold 
M of dimension m and tt : P ^ M be a fibration. We consider the subset of i? x TP 

T^P ^ {{a,v) e E xTP I (r7r)(w) = p{a) } , 

where Tn : TP TM is the tangent map to tt. Denote by r'' : 7^P P the 
map given by T'^{a,v) = Tp{v), Tp : TP — > P being the canonical projection. If 
dimP = p, one may prove that 7^ P is a vector bundle over P of rank n + p — m 
with vector bundle projection t'^ : 7^P P. 

A section X oir^ : T^P ^ P is said to be projectable if there exists a section X 
oi T : E ^ M and a vector field U S X{P) which is 7r-projectable to the vector field 
p{X) and such that X{p) — {X{Tr {p)),U{p)), for all p £ P. For such a projectable 
section X, we will use the following notation X = {X, U). It is easy to prove that 
one may choose a local basis of projectable sections of the space Sec(T'^). 

The vector bundle r'^ : T^P P admits a Lie algebroid structure (|-, ■l'^,p^). 
Indeed, if {X, U) and {Y, V) are projectable sections then 

I(x,c/),(y,i/)r = (ix,y],[c/,F]), p^{x,u) = u. 

(T^P, |-, J^, p'^) is the E-tangent bundle to P or the prolongation of E 
over the fibration tt : P ^ M (for more details, see ^5]). 

Now, let {E, |-, •], p) (resp., (£", |-, •]', p')) be a Lie algebroid over a manifold M 
(resp., M') and suppose that ^ : i? — > i?' is a vector bundle morphism over the map 
^0 '■ M — ^ M' . Then, the pair (^f, ^^q) is said to be a Lie algebroid morphism if 

d((«',^'o)*0') = (^',«'o)*(d' </>'), for all (t>' E Sec(A''XT')*) and for all k, (2.3) 
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where d (resp., d') is the differential of the Lie algebroid E (resp., E') (see [33]). In 
the particular case when M — M' and ^'o = Id then (2.3 1 holds if and only if 

|*oX,*or]' = 5'|X,rl, p'{^X)^p{X), for X,y e Sec(r). 

2.2. Lie groupoids. A Lie groupoid over a differentiable manifold A/ is a differ- 
entiable manifold F together with the following structural maps: 

• A pair of submersions a : T ^ M, the source, and /? : F ^ M, the 
target. The maps a and f3 define the set of composable pairs 

T2 = {{9,h)eGxG \ f3{g) = a{h) } . 

• A multiplication m : F2 ^ F, to be denoted simply by m{g, h) — gh, 
such that 

- aigh) = a{g) and f3{gh) = f3{h). 

- g{hk) = {gh)k. 

• An identity section e : M ^ T such that 

- = 9 and ge{f}{g)) = g. 

• An inversion map i : F ^ F, to be simply denoted by i{g) ~ g^^, such 
that 

- 9^^9 = and gg-^ = e{a{g)). 

A Lie groupoid F over a set M will be simply denoted by the symbol F =^ M. 

On the other hand, if g g F then the left-translation by g and the right- 
translation by g are the diffeomorphisms 

Ig : a^^{f3{g)) — > a-^{a{g)) ; h — > lg{h) ^ gh, 
-^rH/9(ff)) ; h^rg{h)^hg. 

Note that = Ig-i and = rg-i. 

A vector field X on F is said to be left-invariant (resp., right-invariant) if 

it is tangent to the fibers of a (resp., f3) and X{gh) ~ {Thlg){Xh) (resp., X{gh) = 
{Tgrh){X{g))), for {g,h) €T2. 

Now, we will recall the definition of the Lie algebroid associated with F. 

We consider the vector bundle t : Er — > M, whose fiber at a point a; e M is 
{Ey)x = V'e(a;)a = KeriT^^x^a). It is easy to prove that there exists a bijection 
between the space Sec(T) and the set of left-invariant (resp., right-invariant) vector 
fields on F. If X is a section of t : Er M, the corresponding left-invariant (resp., 
right- invariant) vector field on F will be denoted X (resp., X), where 

X{g) = {T,(p(g))lg){X{P{g))), (2.4) 

^(5) = -(T,(„(,))rg)((r,(„(g))z)(X(a(5)))), (2.5) 
for g £ T. Using the above facts, we may introduce a Lie algebroid structure 
(|-, •], p) on Er, which is defined by 



{X,Y\^[X,Y], p{X){x) = {T,^,)p){X{x)), (2.6) 
for X,Y e Sec(T) and x ^M. Note that 



IX,YI = -[X,Y], [X,Y] = Q, (2.7) 

(for more details, see [711^). 

Given two Lie groupoids F =^ M and F' =^ M' , a morphism of Lie groupoids 

is a smooth map <!> : F ^ F' such that 

{g,h)eT2=^mg),<i>{h))e{T% 
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and 

^gh)=^g)^h). 

A morphism of Lie groupoids $ : F F' induces a smooth map : M ^ M' in 
such a way that 

a' o $ = $0 o a, o $ = $0 o /3, <& o e = e' o <I)o, 

a, /3 and e (resp., a', /?' and e') being the source, the target and the identity section 
of F (resp., F'). 

Suppose that ($, $o) is a morphism between the Lie groupoids F =| M and F' =t 
M' and that t : Er M (resp., r' : i5r' ^ M') is the Lie algebroid of F (resp., 
r'). Then, if x € M we may consider the hnear map E^i^) : {Er)x — * {Er')<!>o{x) 
defined by 

Ex{<i>){v^(^x^) = (Te(^)$)(w,(2;)), for e (-Br):.. (2.8) 

In fact, we have that the pair {E{^), <i>o) is a morphism between the Lie algebroids 
t: Er ^ M and r' : Er' M' (see [25]). 

Trivial examples of Lie groupoids are Lie groups and the pair or banal groupoid 
M X M, M being an arbitrary smooth manifold. The Lie algebroid of a Lie group 
r is just the Lie algebra q of F. On the other hand, the Lie algebroid of the pair 
(or banal) groupoid M x M is the tangent bundle TM to M . 

Apart from the Lie algebroid -Er associated with a Lie groupoid F =^ Af , other 
interesting Lie algebroids associated with F are the following ones: 

• The -Er- tangent bundle to Ep: 

Let T^r^;* I-,,;; Er-tangent bundle to Ep, that is, 

T^^E* = { {vx,XrJ e (Er). X Tx^E* | (Tt,t*)(XtJ = (r,(,)/3)(t-,) } 

for Ta, e (Ep)a;, with x e M. As we know, 7^^ E^ is a Lie algebroid over Ep. 

We may introduce the canonical section O of the vector bundle {7^^ E^)* Ep 
as follows: 

for G (-^'r)^ and {ax,Xx^) G T^'^Ep. O is called the Liouville section as- 
sociated with Ep- Moreover, we define the canonical symplectic section 

associated with Er by U — —dQ, where d is the differential on the Lie algebroid 
is easy to prove that is nondegenerate and closed, that is, it is 
a symplectic section of T^"^ Ef (see ^ ) . 

Now, if Z is a section of r : Er M then there is a unique vector field Z*"^ on 
Ep, the complete lift of X to Ep, satisfying the two following conditions: 

(i) Z*'^ is T*-projectable on p{Z) and 

(ii) Z*%X) = (Z^j 

for X G Sec(r) (see |13]). Here, if X is a section of r : Er M then X is the 
linear function X G C°° (E* ) defined by 

X{a*) = a*{X{T*{a*))), for all a* G E*. 

Using the vector field Z*"^, one may introduce the complete lift Z*'^ of Z as the 
section of r""' : 7^^El E^ defined by 

Z*''{a*) = {Z{T*{a*)),Z*''{a*)), for a* G E* . (2.9) 
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Z*'^ is just the Hamiltonian section of Z with respect to the canonical symplectic 
section O associated with Ey. In other words, 

iz.cr^^dZ, (2.10) 

where d is the differential of the Lie algebroid t'^ : Ep (for more details, 

see P]V 

• The Lie algebroid tt : T^F T : 

Let T'^r be the Whitney sum VP (Br Va of the vector bundles VP T and 
Va — > r, where Vf3 (respectively, Va) is the vertical bundle of f3 (respectively, 
a). Then, the vector bundle rp : = V/3 S)r Vet T admits a Lie algebroid 
structure (|-, -J^ ^ , The anchor map p'^ ^ is given by 

and the Lie bracket bracket |-,-]^ ^ on the space Sec(rr) is characterized for the 
following relation 



lix,Y),{rx)f ^^{-lx,x%iY, y '1 ) , 

for X,Y,X',Y' e Sec(r) (for more details, see [27]). 

On other hand, if X is a section of r : Ep M, one may define the sections 
^(1,0)^ jj^(o,i) (the /3 and a-lifts) and X(i'i) (the complete lift) ofXtofr : T^F ^ F 
as follows: 

X(i'°)(5) = (^(5), 0,), X(°'i)(g) = {Og,X{g)), and X^^^^^g) = {-X{g),X{g)). 
We have that 

[x(i'O), r(i'0)f - -fx, ri(i'O) y(i^o)f ^ 

and, as a consequence. 

Now, if g,h Cz T one may introduce the linear monomorphisms ^^'^"^ : (^'r)* 
(T^F)* ^ V^P © y^Ta and : (Sr);^^) ^ (T^r)* ^ V; P ® V^a given by 

Jl'^°\x^„Y^,)=JiThiior^-^){X^,)), (2.11) 

7fi)(^.,i^.)=7((T,Z,-i)(r,)), (2.12) 
for {Xg,Yg) e T^F and (X,„y,0 e T^F. 

Thus, if /i is a section of r* : i?p M, one may define the corresponding lifts 
^(1.0) and fi^-^'^^ as the sections of fr* : (T'^F)* F given by 

p(i.")(/,,) = ^(^i^o\ for/ieF, 
^(0-i)(5) = Mr\ for.9eF. 

Note that if g G F and {Xa} (respectively, {^5}) is a local basis of Sec(T) on an 
open subset U (respectively, V) of M such that a{g) G U (respectively, P{g) G V) 
then Y^°'^^} is a local basis of Sec(rr) on the open subset a~^{U)np~^{V). 

In addition, if {X"^} (respectively, {Y^}) is the dual basis of {Xa} (respectively, 
{Yb}) then {(X'4)(i^o) (yB)(o,i)} the dual basis of yj"'^^}. 
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2.3. Discrete Unconstrained Lagrangian Systems. (See [27] for details) A 
discrete unconstrained Lagrangian system on a Lie groupoid consists of a 
Lie groupoid T =t M (the discrete space) and a discrete Lagrangian L^t :V ^ 
R. 

2.3.1. Discrete unconstrained Euler-Lagrange equations. An admissible sequence 
of order N on the Lie groupoid F is an element (51, ... , qn) of = F x • • • x F 
such that {gk, gk+i) e Fj, for fc = 1, . . . , TV - 1. 

An admissible sequence {gi, . . . ,gN) of order is a solution of the discrete 
unconstrained Euler-Lagrange equations for if 

d''[Ld o Ig^ +LdO rg^^^ o i]{e{xk))\(Er)^^ = 

where P{gk) — <^(5fc+i) — and 0?° is the standard differential on F, that is, the 
differential of the Lie algcbroid Tr : TT — > F (see [27]). 

The discrete unconstrained Euler-Lagrange operator D del Ld ■ F2 Ep 

is given by 

{DDELLd){g,h) = d°[Ld olg + Ldo r,i o i]{e{x))i(Er)^ = 0^ 

for (5, h) £ F2, with f3{g) = a{h) ^xe M (see [27]). 

Thus, an admissible sequence (gi, . . . , gN) of order A^ is a solution of the discrete 
unconstrained Euler-Lagrange equations if and only if 

{DDELLd)igk,gk+i) = 0, for fc = 1, . . . , A^ - 1. 

2.3.2. Discrete Poincare-Cartan sections. Consider the Lie algebroid Tp : T'^F = 
VP (Br Vet F, and define the Poincare-Cartan 1-sections 0^^,Gj^ G 
Sec((-nr)*) as follows 

Ql^{g){Xg,Yg)^-Xg{Ld), e+^{g){Xg,Yg)^Yg{Ld), (2.13) 

for each 5 G F and (Xg^Yg) e T^F = VgP ® Vga. 

Since dLd = Qt — ©7 and so, using = 0, it follows that dQt = dQj . This 
means that there exists a unique 2-section — —dQ^^ = —dQ^^, which will be 
called the Poincare-Cartan 2-section. This 2-section will be important to study 
the symplectic character of the discrete unconstrained Euler-Lagrange equations. 

If g is an element of F such that a{g) = x and (3{g) = y and {A^i} (respectively, 
{Yb}) is a local basis of Sec(r) on the open subset U (respectively, V) of M, with 
X £ U (respectively, y E V), then on a^^{U) n (3^^{V) we have that 

e-Ld = -^a{l){x^Y'''\ el = rEiL)iY^r^), 

^Ld = -Al(l^(Lrf))(A^)(i-")A(r^)(0'i) ^' ' 

where {A'^} (respectively, {Y^}) is the dual basis of {A^} (respectively, {ye}) 
(for more details, see [27]). 

2.3.3. Discrete unconstrained Lagrangian evolution operator. Let T : F — )■ F be a 
smooth map such that: 

- graph(T) C F2, that is, (g, T{g)) £ F2, for all g G F (T is a second order 
operator) and 

- {g,T{g)) is a solution of the discrete unconstrained Euler-Lagrange equa- 
tions, for all g £T, that is, {DoELLd){g, Tr(g)) — 0, for all g £ F. 
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In such a case 

X{g){L,)~X{T{g)){Ld) = 0, (2.15) 

for every section X of r : Er M and every 5 S F. The map T : F ^ F is called 
a discrete flow or a discrete unconstrained Lagrangian evolution operator 
for Ld- 

Now, let T : F ^ F be a second order operator. Then, the prolongation TT : 
T^F = V/3 ®r Va 7^T = V(3 ©r Va of T is the Lie algcbroid niorphism over 
T : F F defined as follows (see [571): 

7J{X,,Y,) = ((r,(r,x(g)Oz))(y,),(T,T)(X3) 

+(^,T)(y,)-T,(r,x(,)o^)(r,)), ^"-^'^ 

for all {Xg,Yg) e T^F = Vg(3 ® Vga. Moreover, from (O, and ( |2.16| , we 

obtain that 

7gT{X{g),Y{9)) = {-Y{T{g)), iTgT){X{g) + Y [g)) + Y{T{g))), (2.17) 
for all X,y G Sec(T). 

Using ( 2.16[ ), one may prove that (see [37]): 

(i) The map T is a discrete unconstrained Lagrangian evolution operator for 
Ld if and only if (TT, T)*e^^ - 9+ . 

(ii) The map T is a discrete unconstrained Lagrangian evolution operator for 
Ld if and only if (TT, T)*e^_^ - = dLd- 

(iii) If T is discrete unconstrained Lagrangian evolution operator then 

(TT,T)*f^i, 

2.3.4. Discrete unconstrained Legendre transformations. Given a Lagrangian Ld '■ 
F ^ M we define the discrete unconstrained Legendre transformations 

F-Ld -.T ^ E* and ¥+Ld -.T E*hy (see [27 ) 

(F"Ld)(/i)(We(a(/i))) = -Ve(a{h}){Ld OThOi), for ^^(^(/i)) G {Er)a{h), 

{¥+ Ld){g){ve(f3(g))) = v^(^fj(^g)){Ld o Ig), for G (£^r)/3(g)- 

Now, we introduce the prolongations T^F^L^ : T^F = VP ©r Va 7^^ E^ and 
T^F+id : T^F = VP ©r Va ^ E^ by 

7l¥-Ld{Xh,Yh) = {n{tor,,-i){Xh),{n¥-Ld){Xh) + {n¥-Ld){n)),{2.18) 

7^g¥ + Ld{Xg,Yg) ={{Tglg-.){Yg),{Tg¥+Ld){Xg) + {Tg¥+Ld){Yg)), (2.19) 

for all h,g er and {Xh,Yh) G T^F = Vhf3 © V^a and (Xg,yg) G T^F = Vg(3 © 
V^a (see [27]). We observe that the discrete Poincare-Cartan 1-sections and 2- 
section are related to the canonical Liouville section of {7^^ E^)* Ep and the 
canonical symplectic section of A^(T^'"£'p)* E^ by pull-back under the discrete 
unconstrained Legendre transformations, that is (see [27]), 

{7^¥- Ld,¥-Ldr Q = Ql^, {7^¥+ Ld,¥+ LdTQ = Q^^, (2.20) 

{7^¥- Ld,¥-Ldr n^n^,, {7^¥+Ld,¥+Ldrn = nL,. (2.21) 

2.3.5. Discrete regular Lagrangians. A discrete Lagrangian : F ^ M is said to be 
regular if the Poincare-Cartan 2-section is nondegenerate on the Lie algebroid 
Tr : T^F = V(3 ©r Va F (see [27]). In [27], we obtained some necessary and 
sufficient conditions for a discrete Lagrangian on a Lie groupoid F to be regular 
that we summarize as follows: 

Ld is regular -4=> The Legendre transformation ¥^ Ld is a local diffeomorphism 
The Legendre transformation ¥^ Ld is a local diffeomorphism 
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Locally, we deduce that L4 is regular if and only if for every 17 e F and every 
local basis {Xa} (respectively, {Yb}) of Sec(T) on an open subset U (respectively, 
V) of M such that a{g) G U (respectively, (3{g) G V) we have that the matrix 
(X a{Y B{.Ld))) is regular on a-\U) n P~\V). 

Now, let : r ^ be a discrete Lagrangian and g be a point of V. We define 
the M-bilincar map G^"* : {Er)a{g) ® (£^r)/3(g) ^ K given by 

G^-{a, b) = r!L,(5)((-r,(„(g))(r, o t){a),0), (0, (T,(^(,));j(6))). (2.22) 



Then, using (2.14), we have that 

Proposition 2.1. The discrete Lagrangian Ld '■ T ^ R is regular if and only if 
Gg'' is nondegenerate, for all g £ T, that is, 

G^''{a,b) = 0, for all b G (i^r)/3(g) ^ a = 

(respectively, Gg'^{a, b) — 0, for all a E {Er)a{g) ^ b = 0). 

On the other hand, if : F ^ M is a discrete Lagrangian on a Lie groupoid F 
then we have that 



where r* : E'^ M is the vector bundle projection. Using these facts, (2.181 and 
(2.191, we deduce the following result. 

Proposition 2.2. Let Ld '■ T —t be a discrete Lagrangian function. Then, the 
following conditions are equivalent: 

(i) Ld is regular. 

(ii) The linear map Ld ■ VhP V^a —> '^r- La(h)-^r linear isomor- 
phism, for all h GT. 

(iii) The linear map 7^¥+Ld : Vg(3 ® VgU -> 7^j:^^^g^Er* is a linear isomor- 



phism, for all g £T. 



Finally, let : F — > M be a regular discrete Lagrangian function and {go, Hq) € 
F X F be a solution of the discrete Euler-Lagrange equations for Ld- Then, one may 
prove (see [23) that there exist two open subsets Uq and Vq of F, with g^ G Uq 
and Hq £ Vq, and there exists a (local) discrete unconstrained Lagrangian evolution 
operator Tl^ : Uq ^ Vq such that: 

(i) TLrf(go) = ^0, 

(ii) T is a diffeomorphism and 

(iii) T^^ is unique, that is, if Uq is an open subset of F, with go G Uq, and 
T'^^ : C/q ^ F is a (local) discrete Lagrangian evolution operator then 

3. Discrete Nonholonomic (or constrained) Lagrangian systems on 

Lie groupoids 

3.1. Discrete Generalized Holder's principle. Let F be a Lie groupoid with 
structural maps 

a,(3 -.T ^ M, e : M ^ F, i:T ^T, m : F2 ^ F. 

Denote by r : Er M the Lie algebroid associated to F. Suppose that the rank 
of Er is n and that the dimension of M is m. 

A generalized discrete nonholonomic (or constrained) Lagrangian system on F 
is determined by: 
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- a regular discrete Lagrangian : T — > K, 

- a constraint distribution, Dc, which is a vector subbundle of the bundle 

Ey- M oi admissible directions. We will denote by tt),, : Tt^ M the 
vector bundle projection and by ij)^ : 1)^ the canonical inclusion. 

- a discrete constraint embedded submanifold Mc of F, such that 
dimMc — dimDc = m + r, with r < n. We will denote by ij^^ : Mc T 
the canonical inclusion. 

Remark 3.1. Let : F ^ M be a regular discrete Lagrangian on a Lie groupoid 
F and Mc be a submanifold of F such that e(M) C Mc- Then, dimMc = m + r, 
with < r < TO. Moreover, for every x G M, we may introduce the subspace Dc{x) 
of Er{x) given by 

Dc(x) = T,(^)Mcni;r(a;). 
Since the linear map T^(x)Ct : T^(^x)Mc T^M is an epimorphism, we deduce that 
dimDc{x) = r. In fact, Dc = {J^eM ^c(a;) is a vector subbundle of Er (over M) of 
rank r. Thus, we may consider the discrete nonholonomic system (Ld,Mc,23c) on 
the Lie groupoid F. o 

For g gT fixed, we consider the following set of admissible sequences of order 

N: 

= { (5fi,...,5Jv) e F^ I (fiffe,5fc+i) e F2, ioT k = l,..,N -1 a,nd gi...gN = g} . 

Given a tangent vector at {gi, . . . ,gN) to the manifold 6^, we may write it as the 
tangent vector at f = of a curve in C^, t e (—£,£) C M — > c{t) which passes 
through ((/I, . . . , giq) at t = 0. This type of curves is of the form 

c{t) = {gihi{t),h^^{t)g2h2{t), h'^\^{t)gN-ihN-i{t),K^\{t)gN) 

where hk{t) G a''^ {I3{gk)) , for all t, and /ifc(O) = e{P{gk)) for = 1, . . . , iV - 1. 
Therefore, we may identify the tangent space to at (51, ... , g^) with 

^(si,S2,..,s«)Gf = {{vi,V2,...,vn-i) I Vk e (£^r)xfc and Xk = /3(5fe),l < fc < A'' - 1 } . 

Observe that each Vk is the tangent vector to the curve hk &t t = 0. 

The curve c is called a variation of (51, . . . , g'jv) and (wi, W2, • • • , ^^iv-i) is called 
an infinitesimal variation of (gi, . . . , gjv). 

Now, we define the discrete action sum associated to the discrete Lagrangian 
Ld-.T — >Ras 



SLd : ef 



s 



N 



{9i,---,9n) I — > ^Li{gk). 



k=l 

We define the variation SSLj, : T(g^^...^g^)e^ M as 
5§Ld{vi,...,VN-i) = ^ §id(c(t)) 



{Ld{gihi{t)) + Ld{h^\t)g2h2{t)) 



t=o 



__ d 
~ dt 

+ . . . + Ld{hJj\^{t)gN-ihN-i{t)) + Ld{hJj^_^{t)gN)} 
= J2 ° l9.Mxk)){vk) + d"{Ld o Tg,^, o i){e{xk)){vk)) 

k=l 

where d° is the standard differential on F, i.e., d° is the differential of the Lie 
algebroid rr : TV — > F. It is obvious from the last expression that the definition 
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of variation dSL^ does not depend on the choice of variations c of the sequence g 
whose infinitesimal variation is {vi,. . . , wjv-i)- 

Next, we will introduce the subset (Vc)g of T^gi,...,gff) 

defined by 

{yc)g = {{vu...,VN-i)eT^g„...,g^)e^ I Vfce {!,..., AT- 1}, VkG-D,}. 

Then, we will say that a sequence in satisfying the constraints determined 
by Mc is a Holder-critical point of the discrete action sum SL^ if the restriction 
of 5$Lci to (Vc)g vanishes, i.e. 



6§Ld 



= 0. 



Definition 3.2 (Discrete Holder's principle). Given g G F. a sequence (gi, .... g.v) 
€ such that g^ G Mg, 1 < k < N, is a solution of the discrete nonholo- 
nomic Lagrangian system determined by {Ld, Mc, ©c) */ md only if {gi,. . . , g^) is 
a Holder- critical point ofSLj^. 

If (51, . . . ,gN) € n (Mc X • • • X Mc) then (51, . . . ,5jv) is a solution of the 
nonholonomic discrete Lagrangian system if and only if 

iV-l 

('^"(^'^ ° ^9.) + d%Ld o rg,^^ o i)){e{xk))iiv^u^ = 0, 

k=l 

where f3{gk) = a{gk+i) = x^. For = 2, we obtain that {g, /i) e r2 fl (Mc x Mc) 
(with P{g) = a{h) = a;) is a solution if 

d°{Ld olg + LdorhO i)(e(x))|(D^)^ = 0. 

These equations will be called the discrete nonholonomic Euler- Lagrange 
equations for the system (L^, Mc, 2)c). 

Let {gi, . . . ,gN) be an element of G^. Suppose that P{gk) = ct{gk+i) = Xk, 
1 < k < N — 1, and that {X^ife} = {Xak, Xak} is a local adapted basis of Scc(t) 
on an open subset Uk of M, with Xk S Uk- Here, {Xak}i<a<r is a local basis of 
Sec(rD^) and, thus, {X"*'}r+i<a<n is a local basis of the space of sections of the 
vector subbundlc tdo : 1)° ^ M, where Dj? is the annihilator of Dc and {X"*", X"*"} 
is the dual basis of {Xak, Xak}- Then, the sequence {gi, . . . ,gN) is a solution of 
the discrete nonholonomic equations if (gi, . . . , gjv) € Mc x • • • x Mc and it satisfies 
the following closed system of difference equations 



= Y ^Xak{gk){Ld) - Xak{gk+l){Ld) , 

k=l 

= J2 [(d^d^ (^afc)("'^))(5fe) - {dLd, (X.fe)(l'0))(5fc+l) 



fe=l 

for 1 < a < r, d being the differential of the Lie algcbroid tt"^ : T'^F = Vf3(BrVa — > 
F. For = 2 we obtain that {g, /i) G F2 n (Mc x Mc) (with P{g) = a{h) = a;) is a 
solution if 

X'a{g){La)~xl{h){Ld) = Q 

where {A^a} is a local basis of Scc(ti)^) on an open subset U of M such that x G U . 

Next, we describe an alternative version of these difference equations. First 
observe that using the Lagrange multipliers the discrete nonholonomic equations 
are rewritten as 

d" [La olg + Laorhoi] (e(a;))(t;) = A„X«(ar)(^;), 
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for V G {Er)x, with {g, /i) e Ta n (Mc x Mc) and j3{g) = a{h) = x. Here, {X"} is 
a local basis of sections of the annihilator D°. 

Thus, the discrete nonholonomic equations are: 

y (ff)(Lrf) - Y{h){La) = A„(x«)(y)|^(,), (5, /i) e r2 n (M^ x mj, 

for all y e Sec(T) or, alternatively, 

(rfirf - A„(X«)(0'i),y(0'i))(ff) - (dLd,y(i'0))(/i) = 0, {g,h) e r2 n (Me x Me), 
for all y e Sec(r). 

On the other hand, wc may define the discrete nonholonomic Euler- Lagrange 
operator DDEL{Ld,Mc,'Dc) : T2 n (Me x Me) D* as follows 

DoELiLd, Me, 2)e)(fl', h) = (F [Ld olg + Ldorhoi] (e(a;)) , 

for (5, /i) e r2 n (Mc x Me), with l3{g) = a{h) = xe M. 

Then, we may characterize the solutions of the discrete nonholonomic equations 
as the sequences ((/i, . . . , gjv), with {gk,gk+i) G Fa fl (Me x Me), for each k G 
{!,..., TV - 1}, and 

DDEL{Ld,'Mc,Dc){9k,gk+i) =0. 
Remeirk 3.3. (i) The set Fa fl (Mc x Me) is not, in general, a submanifold 

of Mc X Me- 

(ii) Suppose that ajvtc '■ M and (3^^ '■ M are the restrictions 

to Mc of a : r ^ Af and f3 : T ^ M, respectively. If ajvtc I^Mc ^-rc 
submersions then T2 fl (Mc x Mc) is a submanifold of Mc x Me of dimension 
m + 2r. 

o 

3.2. Discrete Nonholonomic Legendre treinsformations. Let {Ld, Me, Dc) be 
a discrete nonholonomic Lagrangian system. We define the discrete nonholo- 
nomic Legendre transformations 

F-(Ld,Me,2)c) :Me-^D* and ¥+{Ld,M^,'D^) : ^ 'D* 

as follows: 

¥-{Ld,Mc,T)c){h){v^f^aih))) = -Ve{a{h)){Ld°rhoi), for Ue(a(/i)) e Dc{o:{h)), (3.1) 
F+(Lrf,Mc,2)e)(.9)(^^cai(<,))) =«e(/3(s))(ido;3), for G ©e(/3(<7)) ■ (3-2) 

If V~Ld : r — ^ i?p and F+L^ : F ^ are the standard Legendre transformations 
associated with the Lagrangian function Ld and : D* is the dual map of 

the canonical inclusion : T>c — > -Br then 

F-(id,Me,Dc) =«B,oF-Ldoi3vc^, ¥+{Ld,Mc,'Dc) =i*^^o¥+LdoiM,. (3.3) 

Remeirk 3.4. (i) Note that 

r^,^oF-(Ld,Me,De) = aM„ t^^ o¥+(Ld,Mc,D^) = (3m^. (3.4) 

(ii) If DDEL{Ld,'^c,'^c) is the discrete nonholonomic Euler-Lagrange opera- 
tor then 

DdEL {Ld, Me, Be) (5, /l) = F+ {Ld, Me, De) (5) ' (^d, Me, ©c) {k) , 

for (5,/!) e F2 n (Me X Mc). 

o 
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On the other hand, since by assumption L^; : F — + M is a regular discrete La- 
grangian function, we have that the discrete Poincare-Cartan 2-section H.^^ is sym- 
plectic on the Lie algebroid : 7^T ^ T. Moreover, the regularity of L is equiv- 
alent to the fact that the Legendre transformations F~Lrf and F+L^; to be local 
diffeomorphisms (see Subsection 2.3.51. 

Next, we will obtain necessary and sufficient conditions for the discrete non- 
holonomic Legendre transformations associated with the system {Ld, Mc, 2)c) to be 
local diffeomorphisms. 

Let F be the vector subbundle (over F) of fp : T'^F F whose fiber at the point 
heT is 

In other words, 

^^,° = {7^°^ I 7 e D.(a(/.))° } . 

Note that the rank of F is n + r. 

We also consider the vector subbundle F (over F) of -rr : T'^F ^ F of rank n 
whose fiber at the point g G F is 



CT^F. 



{7^^ I 7 e B,(/3(g))" 



CT^^F. 



Lemma 3.5. F (respectively, F) is a coisotropic vector subbundle oj the symplectic 
vector bundle (T'^Fjili^), that is, 

Ff^ C Fin for every h £T 

(respectively, F^ C Fg, for every g & T), where F^ (respectively, F^) is the 
symplectic orthogonal of F^ (respectively, Fg) in the symplectic vector space (T^j^F, 
^lA^)) (respectively, (T^F, 



Proof. If e F we have that 



TjjF (TjjF)* being the canonical isomorphism induced by the symplectic 



form ilL^{h). Thus, using (2.141, we deduce that 
F, 



h - { Ku^iH)ilH''^) I 7 e B,(a(/.))" } C {0} ® V^a C F,. 
The coisotropic character of Fg is proved in a similar way. □ 

We also have the following result 

Lemma 3.6. Letl^W-Ld : T^F ^ T^^^* (respectively, T^F+id : T^^F ^ 7^^E^) 
be the prolongation of the Legendre transformation ¥^ L^i : F E^ (respectively, 
¥+Ld -.T ^ E^). Then, 



{7l¥-Ld){Fn) = 7, 
for h G Mc (respectively, 

{7l¥+Ld)iFg) = T,^+^^(^)F? = { («;3(,),XF+i,(,)) e 
for g(EMJ. 



Ld(g) i 



Va{h) e T)c{a{h)) 
vpig) e Dc(/3(5)) } 



Proof It follows using ( |2.11| ), ( |2.18[ ) (respectively, ( |2.12[ ), ( |2.19[ )) and Proposition 



□ 



Now, we may prove the following theorem. 



16 D. IGLESIAS, J. C. MARRERO, D. MARTfN DE DIEGO, AND E. MARTfNEZ 



Theorem 3.7. Let (Ld,M.c,Dc) be a discrete nonholonomic Lagrangian system. 
Then, the following conditions are equivalent: 

(i) The discrete nonholonomic Legendre transformation V~ {Ld,'M.c,T)c) (re- 
spectively, F+(Ld, Mc, ©c) j is a local diffeomorphism. 

(ii) For every h € Mc (respectively, g € JAc) 



{p-^'^r^TMnF^^io} 

(respectively, {p'^''^)-^{TgMc) n {0};. 



(3.5) 



Proof (i) => (ii) If heMc and (X,„ Y^) e (p^ ^)~\T,Mc) n F^ then, using the 
fact that Ff-^ C {0} Vha (see the proof of Lemma 3.5 1, we have that X/^ — 0. 
Therefore, 

y,, e Vha n T,,M,. (3.6) 

Next, we will see that 

in¥-iLd,Mc,^c)){Yh)^0. (3.7) 



From (|3^ and (|3^, it follows that {Th¥- {Ld,Mc,1)c)){Yh) is vertical with 
respect to the projection : CD* ^ M. 

Thus, it is sufficient to prove that 

i{n¥-{La,Mc,V,)){Yh)){Z) = 0, for all Z e Sec(Ta,J. 
Here, Z : D* ^ M is the linear function on D* induced by the section Z. 



Now, using (3.3 1, we deduce that 
((r„F-(Lrf,M„2),))(y;,))(Z)=d(ZoiiJ((F-Lrf)(/i))(0,(T^F-Ld)(n)), 

where d is the differential of the Lie algebroid t"^ : 7^'^ Ep E^. 

Consequently, if : Ep E^ is the complete lift oi Z E Sec(T), we have 

that (see (|2.10[)), 



((r^F-(Ld,Mc,2)c))(l^/0)(^) = f^(F-L<i(/i))(Z*'=(F-Ld(/i)), 

(o,(T,F-L,)(ro), 



(3.8) 



being the canonical symplectic section associated with the Lie algebroid Er- 
On the other hand, since Z 6 Sec(rD^), it follows that Z*'^(¥^ Lij{h)) is in 

we conclude that there exists (X^,Y^) £ Fh 



•J T 



¥-La{h:)^r 
such that 



ii'p and, from Lemma 



3.6 



{7l¥-Ld){XlYl) = Z*^{{¥-L,){h)). 
Moreover, using ( 2.18[ ), we obtain that 

{7l¥-Ld){0,Yh) - {0,{n¥-Ld){n)). 



(3.9) 
(3.10) 



Thus, from ( |2.21| , ( |3.8[ ), ( |3.9| and ( |3.10| , we deduce that 

((r,F-(irf,M,2)e))(n))(^) = -r!i,(/j)((o,y;,),(x;,r,:)). 



Therefore, since (0, Yh) G F^, it follows that (3.7) holds, which implies that Y^ = 0. 
This proves that {p'^''^)-^{ThMc) n F,-^ = {0}. 

If F+(L(i, Mc, 2)c) is a local diffeomorphism then, proceeding as above, we have 
that (p^'T)-i(rgMc) n {0}, for all g G M^. 

(ii) ^ (i) Suppose that h G Mc and that is a tangent vector to Mc at h such 
that 

(r;,F-(Ld,Mc,Dc))(lh) = 0. (3.11) 
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We have that {Tha){Yfi) = and, thus, 

(0,n) e {p'^"^)-\nMc). 

± 

h 



We will see that (0, Yh) G F^, that is 



f]L,(/i)((o,y^), {xiYi)) = 0, for (x,;,y^) g (3.12) 

Now, using ( |2.18[ ) and (2.21), we deduce that 
n^Ahm. Y„), (XlYD) - n(¥-LAh))m (T,F-L,)(y,0), {7l¥-L,){Xl Y^)). 
Therefore, from Lemma |3.6| we obtain that 

f]i,(/i)((o,y;,), - f^(F-L,(/i))(o, (r;,F-irf)(n)), Kw,>F-L.(h))) 

with e 7^_^^^^f^^E^. 

Next, we take a section Z E Sec(Ti)^) such that Z{a{h)) = Va(h)- Then (see 



where Y^- La(h) ^ -^f^idCi)-^? ^^"^ ^F-Ld(h) vertical with respect to the projection 
T* -.E^^ M. 

Thus, since (see Eq. (3.7) in [23]) 

f)(F-X,(/i))((0,(T,F-L,)(n)),(0,y;_^^(,))) = 0, 

we have that 

nLAh)m Yu), {XI y^)) = -n{v-LAh)){z*Av-La{h)), (0, (r;,F-i,)(y,0)) 

= -d{Zoi*^ ){¥-La{h)){Q,{Th¥-Ld){Yh)) 



and, from (3.111, we deduce that (3.12) holds. 

This proves that Yu G {p^^^y^iThMc) n which implies that Yh = 0. 
Therefore, F~(L(i, Mc, 2)c) is a local diffeomorphism. 

If {p^''^)~^{Tg'Mc) n Fg"- = {0} for all g € Mc then, proceeding as above, we 
obtain that F+(Ld, Mc, Dc) is a local diffeomorphism. □ 

Now, let p^^^ : T^F TV be the anchor map of the Lie algebroid tt^ : T^F T. 
Then, we will denote by 'Kt the subspace of T^^^F given by 

= (p^^^)-i(T,,Mc) n Fh, for h e Me. 
In a similar way, for every g g Mc we will introduce the subspace 'Kg of F defined 

by 

% = {p'"^)-\Tgl^c)r,Fg. 

On the other hand, let ft, be a point of Mc and G^'' : {Er)a{h) © (£'r)/3(/i) — > K 
be the K-bilinear map given by (2.22). We will denote by (i?r)^'' the subspace of 
iEr)i3(h) defined by 

{Er)f^ ^{be (i?r)/3(h) I (T^imMb) e UMc } 

and by G^"" : (2)c)a(h) x (Er)^" M the restriction to ('Dc)a(h) x (Er)^" of the 
R-bilinear map G^**. 

In a similar way, if 17 is a point of F we will consider the subspace (Er)'^'' of 
{Er)aig) defined by 

(Er)f^ = {ae {Er)o.ig) I iT,^o.(9)){rg o z))(a) € T^M^ } 
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and the restriction G^-"" : {Er)f- x (I'c)/3(g) M of G^" to the space (Sr)f = x 

('Dc)/3(g)- 

Proposition 3.8. Let (L^, Mc, be a discrete nonholonomic Lagrangian system. 
Then, the following conditions are equivalent: 

(i) For every h g Mc (respectively, g G JAc) 

ip^"'')-\nMc)nF,-t^{0} 

(respectively, {p'^''^)-^{TgMa) n = {0};. 

(ii) For every h G Mc (respectively, g G JAc) the dimension of the vector sub- 
space !K/i (respectively, !Kg) is 2r and the restriction to the vector subbundle 
!K (respectively, of the Poincare-Cartan 2-section flj^^ is nondegener- 
ate. 

(iii) For every h £ Mc (respectively, g € JAc) 

{ b e (Er)^^ I G^-^{a, b) - 0, Va e (2),)„(„) } = {0} 

(respectively, { a G (Er)f' | G^-'{a,b) = 0,V5 e (Dc)^(g) } = {0};. 
Proof, (i) (ii) Assume that h E Mc and that 

(p^'i^)-i(r,Mc)nF,^ = {0}. (3.13) 

Let (7 be an open subset of F, with h £ U, and {'/''''}7=i....,n-r a set of independent 
real C°°-functions on U such that 

M^nU ^{h' eU I </>'^(/i') = 0, for all 7 } . 

If d is the differential of the Lie algebroid tt : T'^F F then it is easy to prove 
that 

(/'^)-i(T,Mc) ^<{dr{h)}>' . 

Thus, 

dim{{p'^'''^y\ThMc)) >n + r. (3.14) 



On the other hand, dimF^ = n — r. Therefore, from (3.131 and (3.141, we obtain 
that 

dim{{p'^^^)-\nM,)) = n + r 

and 

7lT^{p'^"^)-\nM,)OF^. 
Consequently, using Lemma |3.5[ we deduce that 

Fh = Kh®F^. (3.15) 



This implies that dimJih — 2r. Moreover, from (3.151, we also get that 

n:Ki<z 'Ku n f^ 

and, since 'Kh n Fj^ = (p'^'^^)"i(r/iMc) n (see Lemma 3.5 1, it follows that 
J{„n?{^ = {0}. 

Thus, we have proved that IK^ is a symplectic subspace of the symplectic vector 
space (JlT.nLAh))- 

If (p^ '")~^(TgMc) n -F^- = {0}, for all g £ Mc then, proceeding as above, we ob- 
tain that 5{g is a symplectic subspace of the symplectic vector space (T^F, ^1^^ (g)), 
for all g e Mc. 

(ii) => (i) Suppose that h £ Mc and that 'Kh is a symplectic subspace of the 
symplectic vector space {7^T,QLdW)- 
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we deduce that 



If {Xh,Yh) e (p^ ^)-^{nMc) n then, using Lemma 

{Xh,Yh)e:Hh. 

Now, if {X'i^,Y/^) E "Kh then, since {Xh,Yh) E Fj^, we conclude that 

f}i,(/i)((x,,r,),(x;,r^)) = o. 

This impHes that 

{Xh,Yh)E^hn:Ki = {0}. 

Therefore, we have proved that {p'^''^)-\Ti,Mc) n Ffj^ = {0}. 

If Jig n JCj- — {0}, for all g G Mc then, proceeding as above, we obtain that 
(p'^^^r^TgMc) n = {0}, for all g G M,. 

(i) (iii) Assume that 

and that b £ (£'r)f^^'' satisfies the following condition 

G^^(a,6) = 0, Vae 
Then, y„ = (r,(^(^))/;,)(fe) e r„M, n ^;,a and (0,r^) e (p^'r)-i(j.^3ytj^ 
Moreover, if e F/j, we have that 

X'h = -{Te(a{h)){rh o i)){a), with a G ('Dc)Q(h)- 



Thus, using ( |2.14[ ) and ( |2.22| , we deduce that 

f^i,(/i)((x;,r;:),(o,y^)) - f^i,(/i)((x;,,o),(o,y,0) = G^(a,6) = o. 

Therefore, 

which implies that b — 0. 

If (p'^'^^y-iTgMc) n F^ ^ {0}, for all g e Mc then, proceeding as above, we 
obtain that 

{ a e (Fr)f ^ I G^-^a, b) = 0, for all b e } - {0}. 

(iii) =^ (i) Suppose that h E Mc, that 

{ b e (Fr)^ I G^^(a, b) = 0, Va e (Dc)c.(/.) } = {0} 

and let {Xh, Yh) be an element of the set {p'^^^)-'^{ThMc) n F,^. 

Then (see the proof of Lemma 3.51, Xh = and Y^ G T^McHV/ia. Consequently, 

Yh = (T,(/3(,,)/^)(6), with e (Fr)f ^ 
Now, if a G ('Dc)q(/i)5 we have that 

K = {Te(o.(h)){rh o i)){a) e and (X;,, 0) G F^. 



Thus, from ( |2.22| and since (0, r,j) G F^, it follows that 
Therefore, b — Q which implies that Y/i = 0. 



If { a e (Fr)f - I G^''=(a,6) = 0,V6 G (Bc)/3(3) } = ^r all .g e Mc, then 
proceeding as above we obtain that (p^ '")^^(TgMc)nF^ = {0}, for all g G Mc. □ 



Using Theorem 3.7 and Proposition |3.8| we conclude 



20 D. IGLESIAS, J. C. MARRERO, D. MARTfN DE DIEGO, AND E. MARTfNEZ 



Theorem 3.9. Let (Ld,M.c,Dc) be a discrete nonholonomic Lagrangian system. 
Then, the following conditions are equivalent: 

(i) The discrete nonholonomic Legendre transformation ¥^ (L^jM-cDc) (re- 
spectively, {L^jJActDc)) is a local difjeomorphism. 

(ii) For every h G Mc (respectively, g G Mc) 

{p^'^)-\TMnFt^{0} 

(respectively, {p'^''^)~^{TgMc) n = {0}). 

(iii) For every h g Mc (respectively, g G "Mc) the dimension of the vector sub- 
space 3-Ch (respectively, !Hg) is 2r and the restriction to the vector subbundle 
!K (respectively, !H) of the Poincare-Cartan 2-section fi^^ is nondegener- 
ate. 

(iv) For every h G Mc (respectively, g G Mcj 

{ b G (Er)^^ I G^^^ia, b) = 0, Va G (2)c)c(„) } - {0} 

(respectively, { a G (Er)f^ | G^--{a,b) = 0,V6 G (Dc)/3(s) } = {0}^ 

3.3. Nonholonomic evolution operators and regular discrete nonholo- 
nomic Lagrangian systems. First of all, we will introduce the definition of a 
nonholonomic evolution operator. 

Definition 3.10. Let (L^jMcBc) be a discrete nonholonomic Lagrangian system 
and Tnh '■ Mc — > Mc be a differentiable map. Tnh is said to be a discrete nonholo- 
nomic evolution operator for (L^jMcDc) if- 

(i) graph{Tnh) ^ '^2, that is, {g,T,^h{g)) G Ta, for all g G Mc and 

(ii) (g,Tnhig)) is a solution of the discrete nonholonomic equations, for all 
g G Mc, that is, 

d°{Ld olg + LdO rT„^(g) o i){e{P{g)))\D^(p(g)) = 0, for all g G Mc. 



Remark 3.11. If Tnh '■ Mc Mc is a differentiable map then, from (3.11, (3.2) 



and ( |3.4| , we deduce that Tnh is a discrete nonholonomic evolution operator for 
{Ld, Mc, T>c) if and only if 

¥-iLd,Mc,Dc)oTnh=V+{Ld,Mc,Dc). 

o 

Now, we will introduce the notion of a regular discrete nonholonomic Lagrangian 
system. 

Definition 3.12. A discrete nonholonomic Lagrangian system (Ld,Mc,2)c) is said 
to be regular if the discrete nonholonomic Legendre transformations F~(Ld,Mc,Dc) 
and F+(_L£;, Mc, ©c) a'^e local diffeomorphims. 

From Theorem 13.91 we deduce 



Corollary 3.13. Let (L^;, Mc, ©c) be a discrete nonholonomic Lagrangian system. 
Then, the following conditions are equivalent: 

(i) The system (L^jMcDc) is regular. 

(ii) The following relations hold 

{p'^^''y\ThM,) n F^ = {0}, for all h G Mc, 
(/^^)-i(rgMc) n F^ = {0}, for all g G Mc- 
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(iii) !K and "K are symplectic subbundles of rank 2r of the symplectic vector 
bundle (T^T, ^lJ- 

(iv) If g and h are points o/Mc then the R-bilinear maps G^"" and G^'^" are 
right and left nondegenerate, respectively. 

The map G^'''^ (respectively, G^'^'^) is right nondegenerate (respectively, left non- 
degenerate) if 

G^''=(a,6)=0,Vae(Dc)oW=^6-0 

(respectively, G^-^^a^b) = 0,V6 e {1^c)p{g} ^ a = 0). 

Every solution of the discrete nonholonomic equations for a regular discrete 
nonholonomic Lagrangian system determines a unique local discrete nonholonomic 
evolution operator. More precisely, we may prove the following result: 

Theorem 3.14. Let (LfjjMcX'c) be o, regular discrete nonholonomic Lagrangian 
system and (gQ, Hq) € Mc x Mc be a solution of the discrete nonholonomic equations 
for (L^;, Mc, Dc)- Then, there exist two open subsets Uq and Vq ofT, with go G Uq 
and ho G Vo, and there exists a local discrete nonholonomic evolution operator 
^1^"^'='^'^ -.UaDMc^VonMc such that: 

(i) Ti^f^-l^Hao) - ho; 
(^^) T^^'*'''^'"^'''' is a difjeomorphism and 

(iii) Y^^'*'''^'"^'''' is unique, that is, if U'o is an open subset ofT, with go € J/g, 
and Tnh ■ U'o C\ Mc Mc is a (local) discrete nonholonomic evolution 
operator then 

l-'-n/i Jiyona^nMc — Un/iJltfoca^nMc- 



Proof. From remark [3. 4[ we deduce that 

F+(L<j,Mc,Dc)(go) =F-(Ld,Mc,Dc)(/io) - Mo e 

Thus, we can choose two open subsets C/q and Vq of F, with go G Uq and /iq G Vq, 
and an open subset Wo of such that no G Wo and 

F+(L<j,Mc,Dc) :C/onMc^ Won©:, F-(Ld, Mc, ©c) : ^o n Mc ^ l^o n 

are diffeomorphisms. Therefore, from Remark |3. 11 we deduce that 

Y(i.„M.,D.) ^ (F-(L<i,Mc,2)c)"'oF+(Ld,Mc,Bc))|c/onM. : C/q n Mc ^ l^o n Mc 

is a (local) discrete nonholonomic evolution operator. Moreover, it is clear that 
^']^h'^''''^'\9n) — h-a and it follows that fi^^'^"''^"^ jg diffeomorphism. 

Finally, if Uq is an open subset of F, with go G Uq, and Tnh ■ U'o fl Mc Mc 
is another (local) discrete nonholonomic evolution operator then {T nh)\Uor\u;jr\Mc 
is also a (local) discrete nonholonomic evolution operator. Consequently, from 
Remark |3.11[ we conclude that 

= [F-(Lrf,Mc,Dc)-ioF+(L<j,Mc,Dc)] 

- y^nh j|(7onC/^nMc- 

□ 
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3.4. Reversible discrete nonholonomic Lagrangian systems. Let {Ld,M.c, 

Dc) be a discrete nonholonomic Lagrangian system on a Lie groupoid F =| M. 

Following the terminology used in j36j for the particular case when F is the pair 
groupoid M x M, we will introduce the following definition 

Definition 3.15. The discrete nonholonomic Lagrangian system (L(i,Mc,2)c) is 
said to be reversible if 

Ldoi^Ld, i{Mc)^Mc, 
i : F — > F being the inversion of the Lie groupoid F. 

For a reversible discrete nonholonomic Lagrangian system we have the following 
result: 

Proposition 3.16. Let (L^i, Mc, ©c) « reversible nonholonomic Lagrangian sys- 
tem on a Lie groupoid F. Then, the following conditions are equivalent: 

(i) The discrete nonholonomic Legendre transformation F^(Lrf, Mc, Dc) is a 
local diffeomorphism. 

(ii) The discrete nonholonomic Legendre transformation F+(_Ld, Mc, Dc) is a 
local diffeomorphism. 



Proof. If /i e Mc then, using (3.1 1 and the fact that LdO i = Ld, it follows that 

¥~{Ld, Mc, 'Dc)(/l)(We(a(/i))) = ~V4a{h)){Ld ° Ih^) 

for W£(Q.(/i)) G (Dc)c((/i)- Thus, from (3.2), we obtain that 

F-(Ld,Mc,Dc)(/l)(fe(a(/.))) = -V+{LdMc,'Dc){h-^){.v,^p^h-r))). 
This implies that 

F+(irf,Mc,2)c) --F-(Lrf,Mc,Dc)oz. 

Therefore, since the inversion is a diffeomorphism (in fact, we have that ~ id), 
we deduce the result □ 

Using Theorem |3.9| Definition |3 . 1 2 1 and Proposition |3.16| we prove the following 
corollaries. 

Corollary 3.17. Let (L^, Mc, Dc) be- o, reversible nonholonomic Lagrangian system 
on a Lie groupoid F. Then, the following conditions are equivalent: 

(i) The system (Ld,M.c,Dc) is regular. 

(ii) For all h e Mc, 

(iii) :K = {p^^)~^{TMc)nF is a symplectic subbundle of the symplectic vector 
bundle (T^^F, ^lJ- 

(iv) The R-bilinear map G^'^'^ : (-Er)^" x {'^c)a{h) — ^ K is right nondegenerate, 
for all /i e Mc. 

Corollary 3.18. Let {Ld, Mc, Dc) be a reversible nonholonomic Lagrangian system 
on a Lie groupoid F. Then, the following conditions are equivalent: 

(i) The system (Ld,Mc,Dc) is regular. 

(ii) For all g G Mc, 

(/'r)-i(r,Mc)n^/ = {0}. 

(iii) ■k = {p'^ ^)-^{TMc)r\F is a symplectic subbundle of the symplectic vector 
bundle (T^^F, ^Lo)- 
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(iv) The R-bilinear map G^"" : (I'c)/3(g) X (Sr)f ^ M is left nondegenerate, 
for all g G Mc- 

Next, we will prove that a reversible nonholonomic Lagrangian system is dynam- 
ically reversible. 

Proposition 3.19. Let (L^jMcDc) be a reversible nonholonomic Lagrangian sys- 
tem on a Lie groupoidT and {g, h) he a solution of the discrete nonholonomic Euler- 
Lagrange equations for (Ld,M.c,1)c)- Then, {h~^,g~^) is also a solution of these 
equations. Ln particular, if the system {Ld, Mc, ©c) is regular and y^^'''''^'='^=^ ig the 

(local) discrete nonholonomic evolution operator for (Ld,Mc,2)c) then T^^/^'''^'"'^'''^ 
is reversible, that is, 

nh nh 

Proof. Using that i(Mc) = Mc, we deduce that 

Now, suppose that (3{g) = a{h) — x and that v e (I'c)x- Then, since Ld o i — Ld, 
it follows that 

d°[Ldolf^-i + LdOTg-i oi]{£{x)){v) = v{LdO io rho i) + v{LdO i o Ig) 

= v{Ld o Ig) + v{Ld o r/i o z) 0. 

Thus, we conclude that {h~^,g~^) is a solution of the discrete nonholonomic 
Euler-Lagrange equations for {Ld,M,c,Dc). 

If the system {Ld, Mc, T>c) is regular and g S Mc, we have that {g, T^^'''^'^''-'(5)) 
is a solution of the discrete nonholonomic Euler-Lagrange equations for {Ld, Mc, Dc)- 
Therefore, (i(T^^'''^'^°''((7)), is also a solution of the dynamical equations 
which implies that 

r^i-'^^''^\^{Tii-^^''^\g)))=^{g). 

□ 



Remark 3.20. Proposition 3.19 was proved in [SBj for the particular case when F 
is the pair groupoid. o 

3.5. Lie groupoid morphisms and reduction. Let (<&, $0) be a Lie groupoid 
morphism between the Lie groupoids F =| M and F' =^ M' . 

Denote by (£^(<I>), <I>o) the corresponding morphism between the Lie algebroids 
Er and Er' of F and F', respectively (see Section 2.2 1. 

If : F ^ M and L'^ : T' ^ M. are discrete Lagrangians on F and F' such that 

Ld = L'dO^ 

then, using Theorem 4.6 in [27], we have that 

{DDELLd){g,h){v) = {DDELL'd){Hg),Hh)){E,{<^){v)) 

for {g, h) e F2 and v G {Et)x, where x = (3{g) = a{h) G M. 
Using this fact, we deduce the following result: 

Corollary 3.21. Let (<i>, <i>o) he a Lie groupoid morphism, hetween the Lie groupoids 
F =^ A/ and F' =^ M' . Suppose that : F' — > M is a discrete Lagrangian on V , 
that {Ld = o <!>, Mc, Dc) is a discrete nonholonomic Lagrangian system on F and 
that {g, ft.) G F2 n (Mc x Mc). Then: 
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(i) The pair [g, h) is a solution of the discrete nonholonomic problem (Ld,'Mc, 
Dc) if and only if (D£iELL[i){^{g),^{h)) vanishes over the set 

(i^/3(,)'i>)((2)c)/3(s)). 

(ii) // (LJ;, M^, 2)^) is a discrete nonholonomic Lagrangian system on V such 
that m9)Mh)) e X and (ii;;3(g)($))((D,)^(,)) = (2)D*o(/3(s)) then 
{g,h) is a solution for the discrete nonholonomic problem (L^, Mc, ©c) 
and only if (<^{g), '^{h)) is a solution for the discrete nonholonomic problem 

3.6. Discrete nonholonomic Hamiltonian evolution operator. Let (Ld,Mc, 
©c) £^ regular discrete nonholonomic system. Assume, without the loss of gener- 
ality, that the discrete nonholonomic Legendre transformations F^(Lci, 3Vtc, 'Dc) : 
Mc — > D* and F+(Ld, Mc, 2)c) : Mc — > D* are global diffeomorphisms. Then, 

Tit"'^"^''' = F~(Ld,Mc,lic)"^oF+(Ld,Mc,Dc) is the discrete nonholonomic evo- 
lution operator and one may define the discrete nonholonomic Hamiltonian 
evolution operator, jnh : D* ^ D*, by 

7„„ =F+(irf,Me,2)c)o7it"^^^^^'°F+(Ld,Mc,2)c)-' • (3.16) 

From Remark |3.11| we have the following alternative definitions 

7„„ = F-(Ld,Mc,Dc)o7^^'''^-^^)oF-(irf,Me,2)e)-\ 
7„,, = F+(Ld,Me,2)c)oF-(L<i,Mc,Dc)"' 

of the discrete Hamiltonian evolution operator. The following commutative diagram 
illustrates the situation 




Remark 3.22. The discrete nonholonomic evolution operator is an application 

from D* to itself. It is remarkable that D* is also the appropriate nonholonomic 
momentum space for a continuous nonholonomic system defined by a Lagrangian 
L : Ey M and the constraint distribution Dc- Therefore, in the regular case, the 
solution of the continuous nonholonomic Lagrangian system also determines a flow 
from D* to itself. We consider that this would be a good starting point to compare 
the discrete and continuous dynamics and eventually to establish a backward error 
analysis for nonholonomic systems. o 

3.7. The discrete nonholonomic momentum map. Let {Ld, Mc, T>c) be a reg- 
ular discrete nonholonomic Lagrangian system on a Lie groupoid F =^ M and 
T : Er ^ M he the Lie algebroid of F. 

Suppose that g is a Lie algebra and that ^l* : g — > Sec(r) is a M-linear map. Then, 
for each x € M, we consider the vector subspace of g given by 



0" = {eefl I 
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and the disjoint union of these vector spaces 



= U 0^ 



We will denote by {q^")* the disjoint union of the dual spaces, that is, 



xeM 



Next, we define the discrete nonholonomic momentum map J"'' : F 
(0^<=)* as follows: J"''(g) G {q^^^'^)* and 



- e+ (vl/(0(i^i))(5) = ^m9){Ld). for g e T and f e g'^^^). 

If ^ : M ^ g is a smooth map such that ^(x) g g^, for all x G AI, then we may 
consider the smooth function : F M defined by 

jnk^g^ = j"''(g)(e(/3(5))), V.g e r. 

Definition 3.23. The Lagrangian Ld is said to be g-invariant with respect ^' if 

^{0^'''HLd) = M^){Ld) - miLd) = 0, ve e g. 

Now, we will prove the following result 

Theorem 3.24. Let x^^'*'"'^'''^''-' : Mc Mc be the local discrete nonholonomic 
evolution operator for the regular system (L^;, Mc, Dc)- If is Q-invariant with 
respect to ^' : g ^ Sec(r) and ^ : M q is a smooth map such that ^(x) e g^, for 
all X G M, then 



for g g Mc. 

Proof. Using that the Lagrangian is g-invariant with respect to we have that 



, (3.17) 

Also, since {g, ^l^h'''''^'"'^''\g)) is a solution of the discrete nonholonomic equations: 



Thus, from (3.171 and (3.181, we find that 

^(l(/3(g))(.9)(L,) = nm9mril'-''-^'\9)){Ld). 
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Therefore, 



[mg))){g){L,) 



□ 

Theorem |3.24| suggests us to introduce the following definition 

Definition 3.25. An element £^ £ g is said to be a horizontal symmetry for 

the discrete nonholonomic system {Ld,JAc, Dc) o,nd the map ^' : g — s- Sec(T) if 

e {T)c)x, for all x G M. 
Now, from Theorem |3.24[ we conclude that 

Corollary 3.26. // is g-invariant with respect to ^' and ^ €z g is a horizontal 
symmetry for (Lfj, Mc, ©c) o,nd ^' : g — > Sec(T) then J^'* : F ^ M is a constant of 
the motion for '^(^^'^"''^"'l ^ ^/j^j^ 

jnh ^ -^(Ld.McDc) _ jnh 

4. Examples 

4.1. Discrete holonomic Lagrangian systems on a Lie groupoid. Let us 

examine the case when the system is subjected to holonomic constraints. 

Let L(i : r — > M be a discrete Lagrangian on a Lie groupoid F =t M. Suppose 
that Mc C F is a Lie subgroupoid of F over M' C M, that is, Mc is a Lie groupoid 
over M' with structural maps 

a|M, : Mc -> Af, : ^ M', e|M' : A^' -> M^, i|M. : -> Mc, 

the canonical inclusions zjvt^ : Mc — » F and iM' '■ M' — > M are injective immer- 
sions and the pair {ij^^, ^m') is a Lie groupoid morphism. We may assume, without 
the loss of generality, that M' = M (in other case, we will replace the Lie groupoid 
F by the Lie subgroupoid F' over M' defined by F' = a-'^{M') n /3-^(Af')). 

Then, if Ljvtc = -^d ° ^'^d rjvt^ : -Emc -'^'-'^ is the Lie algebroid of Mc, we 
have that the discrete (unconstrained) Euler-Lagrange equations for the Lagrangian 
function Ljvtc ^I's: 

X(g)(iMj-^(/i)(iMj-0, (5,/l) e (Me)2, (4.1) 

for X e Sec(T]v[^). 

We are interested in writing these equations in terms of the Lagrangian Ld 
defined on the Lie groupoid F. From Corollary 4.7 (iii) in [27], we deduce that 



{g, h) G (Mc)2 is a solution of Equations 4.1 if and only if DoEhLdig, h) vanishes 



over Im(i5^(g)(i]vt^))- Here, E{ij^J : Ej^^ Er is the Lie algebroid morphism 
induced between Ej^^ and Er by the Lie groupoid morphism (ijvt^ , Id) . Therefore, 
we may consider the discrete holonomic system as the discrete nonholonomic system 

(Ld,Mc,Dc), where V, = (S(%J)(^mJ = Em^- 
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In the particular case, when the subgroupoid Mc is determined by the vanishing 
set oi n — r independent real C°°-functions 0''' : F ^ M: 

Mc = {geT \ 0^(<?) = O, for all 7}, 

then the discrete holonomic equations are equivalent to: 

Y{g){Ld)-Y{h){Ld) = A^d>^(e(/?(g)))(y(/3(ff)), 

= = 0, 

for all Y £ Sec(r), where d° is the standard differential on T. This algorithm is 
a generalization of the Shake algorithm for holonomic systems (see [101 HOI IMl 
for similar results on the pair groupoid Q x Q). 

4.2. Discrete nonholonomic Lagrangian systems on the pair groupoid. 

Let {Ld, Mc, Dc) be a discrete nonholonomic Lagrangian system on the pair group- 
oid Q X Q ^ Q and suppose that (go, gi) is a point of Mc. Then, using the results 



of Section 3.1 we deduce that ((go, 91), (<Zi, 52)) € (Q x Q)2 is a solution of the 



discrete nonholonomic Euler-Lagrange equations for (Ld, Mc, ©c) if and only if 



{D2Ldiqo, qi) + DiLd{qi,q2))\D^(q^) = 0, 
(91,92) e Mc, 



or, equivalently. 



D2Ld{qa,qi) + DiLd{qi,q2) = ^ XjA^qi), 
(91,92) e Mc, 

where Xj are the Lagrange multipliers and {A^} is a local basis of the annihilator 
B". These equations were considered in [10 and [36 . 

Note that if (gi, 92) G L = Q x Q then, in this particular case, G'^^^ : Tq-^Q x 
Tq^Q ^ M is just the M-bilinear map (_D2-Diid)(9i, 92)- 
On the other hand, if (gi, 92) G Mc we have that 



= { e Tq,Q I (0, vq,) e r(,,,,,)Mc } , 



Thus, the system {Ld,'Mc,1>c) is regular if and only if for every (91,92) € Mc 
the following conditions hold: 



If«,, e(TQ)^;^^^)and 

{D2DiLd{qi,q2)Vq^,Vq^) = 0, Vw,2 G Dc(92) 

and 



If«,, e(TQ)^;^,^)and 

{D2DiLd{qi,q2)vq^,Vq.^) = 0, Vw,! G 'Dc(9i) 

The first condition was obtained in in order to guarantee the existence 
of a unique local nonholonomic evolution operator "f ^^^'^'''^'"^''^ for the system 

(L(i,Mc,Dc)- However, in order to assure that x^^^'*'"'^^'^'''' is a (local) diffeomor- 
phism one must assume that the second condition also holds. 
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Example 4.1 {Discrete Nonholonomically Constrained particle). Consider 
the discrete nonholonomic system determined by: 



a) A discrete Lagrangian L, 
, 1 

Ld{xo,yo,zo,xi,yi,zi) = - 



Xi- xo\ , ( yi~yo\ . ( zi- zq-'^ 



h J \ h 

b) A constraint distribution of Q = M^, 

( d d d 

Dc = span iXi^ — + 2/^, ^2 = ^ 
|_ ox oz ay 

c) A discrete constraint submanifold Mc of x determined by the con- 
straint 

X ^1-2^0 ( y\ + yQ\ ( 

0(xo, 2/0,^0,2:1,2/1, 2^1) = — I — ^ — I I — — 

The system (i^, Mc, Dc) is a discretization of a classical continuous nonholonomic 
system: the nonholonomic free particle (for a discussion on this continuous system 
see, for instance, ^1 13)- Note that if i?(H3xB3) — is the Lie algebroid of the 
pair groupoid x =^ then 

T'(a:i,ai,2i,a;i,yi,2i)^c H -E(R3 xB^) (^^i , 2/i , ^i) = I'c (a;i , 2/l , ^1 ) ■ 

Since 

<— 9 d ^ _ ^ ^ 

— ^ 1" 2/17^ — , Ai — — 2/07^ — , A2 — — , A2 — --^ — , 

ax\ az\ oxq ozq ayi oyo 

then, the discrete nonholonomic equations are: 

2/1 ( r2 ) = 0' (4-2) 

0, (4.3) 



2/2 - 22/1 + 2/0 



/l2 

which together with the constraint equation determine a well posed system of dif- 
ference equations. 

We have that 

D2DiLd = —j^{dxo A dxi + dyo A dyi + dzo A dzi} 

(™^)(^o,yo,^o,2:i,ai,zi) = {'^Ogf^ + ^Oal; + Co^ e T(^xo,vo,zo)'^^ I 

Co = 5(00(2/1 + 2/0) - ^0(2:1 - Xq))}. 

ci = 5(01(2/1 + 2/0) + ^1(2:^1 - a;o))}- 
Thus, if we consider the open subset of Mc defined by 

{ (a;o,2/o, ^o,a;i,2/i,zi) G Mc | 2 + ty^ + 2/12/0 7^ 0, 2 + 2/0 + 2/02/1 7^ } 

then in this subset the discrete nonholonomic system is regular. 

Let : = M2 — , x(r3) given by *(a,6) = + 6^. Then g^- = 
span{^'(0 = Ai}, where f : ^ 1^2 jg defined by l(x,y,z) = (1,2/). More- 
over, the Lagrangian is g-invariant with respect to ^ . Therefore, 

^?''(a;i, 2/1,^1, 2^2, 2/2, 2:2) - J|'*(a::o,2/o,zo,a;i,2/i,zi) 

< 

= *(0,2/2 - 2/i)(a:i, 2/1,^1, 2^2, 2/2, 2:2)(id), 
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that is, 



X2~'Xi Z2- Zi\ Xi- Xo Zi- Zo\ . Z2- Zi 



This equation is precisely Equation (4.2 1 



4.3. Discrete nonholonomic Lagrangian systems on a Lie group. Let G be 

a Lie group. G is a Lie groupoid over a single point and the Lie algebra g of G is 
just the Lie algebroid associated with G. 

If g, h £ G, Vh G TfiG and £ T^G we will use the following notation: 

gvh = {Thlg){vh) e TghG, Vhg = {Thrg){vh) e ThgG, 
gat = (T;\/g-i)(aft) e T;^G, ang = {T;^grg-i){ah) &T;^gG. 

Now, let Mc, ©c) be a discrete nonholonomic Lagrangian system on the Lie 
group G, that is, L^; : G ^ M is a discrete Lagrangian, Mc is a submanifold of G 
and Dc is a vector subspace of g. 

If 5i G Mc then (.91,52) G G x G is a solution of the discrete nonholonomic 
Euler-Lagrange equations for {Ld, Mc, Dc) if and only if 

n— r 

9i^dLd{gi) - dLd{g2)g2^ = J2^3^^^' (44) 
5fc G Mc, fc = 1,2 

where Xj are the Lagrange multipliers and {/i-' } is a basis of the annihilator of 
Dc- These equations were obtained in 0^ (see Theorem 3 in [55]). 

Taking = dLd{gk)9k^, fc = 1, 2 then 

n— r 

5fc G Mc,fc = 1,2 

where Ad : G x g — > g is the adjoint action of G on g. These equations were 
obtained in 14J and called discrete Euler-Poincare-Suslov equations. 

On the other hand, from ( 2.14[ ), we have that 



n^Ai^, V), (If', 7r')) = t'(V(id)) - ^{y{Ld))- 

Thus, if g G G then, using ( 2.22 1 , it follows that the M-bihnear map G^ : g x g ^ M 
is given by 

Gg^''(C,^)--V(5)(?(id))- 

Therefore, the system (irf,Mc,2)c) is regular if and only if for every g G Mc the 
following conditions hold: 

V e g/V(5) e T<,Mc and V(5)(T(irf)) = 0,Ve G Dc =^ - 0, 



e G g/ e (.9) G TgMc and V(.9)( ^ (id)) = 0, Vr; G Dc =^ ^ = 0. 
We illustrate this situation with two simple examples previously considered in 
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4.3.1. The discrete Suslov system. (See |14j ) The Suslov system studies the motion 
of a rigid body suspended at its centre of mass under the action of the following 
nonholonomic constraint: the body angular velocity is orthogonal to some fixed 
direction. 

The configuration space is G = SO{3) and the elements of the Lie algebra so(3) 
may be identified with and represented by coordinates (w^jWj^,^^). Without 
loss of generality, let us choose as fixed direction the third vector of the body frame 
ei, 62, 63. Then, the nonholonomic constraint is lOz — 0. 

The discretization of this system is modelled by considering the discrete La- 
grangian Ld : SO{3) — > R defined by Ldi^) = ^Tr (nj), where J represents the 
mass matrix (a symmetric positive-definite matrix with components ( Jij)i<ij<3). 

The constraint submanifold Mc is determined by the constraint Tr (fiSs) = 
(see [Tl]) where 

\ / 1 

^1 = I -1 , = I , ^3 = 

10/ \ -1 

is the standard basis of so(3), the Lie algebra of SO{3). 

The vector subspace T),, — spanjiJi, £^2}- Therefore, — span{E^}. Moreover, 
the exponential map of 5*0(3) is a diffeomorphism from an open subset of "Dc (which 
contains the zero vector) to an open subset of Mc (which contains the identity 
element I). In particular, T/Mc = Dc- 

On the other hand, the discrete Euler-Poincare- Suslov equations are given 

by 

l(r!i)(L<j)-K(r!2)(i<i) = 0, TTin,E3) = 0, ze{l,2}. 

After some straightforward operations, we deduce that the above equations are 
equivalent to: 

Tr {{E,n2 - niE,)J) = 0, Tr {n^E^) = 0, z G {1,2} 
or, considering the components il^ = i^ij'') of the elements of S0{'3), we have that: 

(J23^33 — J33fl'^2 + ^^22^^23^ \ _ f -~J23^'33 — ^^22^^32' ~ '^12^^31^ j 
^^23^22'' + "^12^^13'' — '^13^12'' / \ +'^33^23'' + '^23^22'' + '^13^21'' / 

C—Jl^fl'^^ + Jss^'^i — Ji2^'"23 \ I Jia^'ss + Jl2^'32 + >^11^^3l' \ 

+ ^23^^21'' ~ '^11^13'' + •^13^^11'' J \ ~J33^13 — ^23^12 ~ '^13^^11'' / 

^12'' ~ ^21^' ^12'' ~ ^21'' • 

Moreover, since the discrete Lagrangian verifies that 

Lain) - ^Tr (nj) = ^Tr (nKj) = La{n-') 

and also the constraint satisfies Tr {^lE^) — — Tr (il^^Es), then this discretization 
of the Suslov system is reversible. The regularity condition in f2 e 50(3) is in this 
particular case: 

77 e 5o(3) /Tr {Ei^-qJ) = 0, Tr {E2nf]J) = and Tr {n-rjEs) = =^ 7? = 
It is easy to show that the system is regular in a neighborhood of the identity /. 



DISCRETE NONHOLONOMIC MECHANICS 



31 



4.3.2. The discrete Chaplygin sleigh. (See [12 |T3]) The Chaplygin sleigh system 
describes the motion of a rigid body sliding on a horizontal plane. The body is 
supported at three points, two of which slide freely without friction while the third 
is a knife edge, a constraint that allows no motion orthogonal to this edge (see [H]). 

The configuration space of this system is the group SE{2) of Euclidean motions 
of M? . An element G SE{2) is represented by a matrix 



cos t) 


— sm 


sin 9 


cos 6 









f2 = I sin cos 9 y | with 9,x,y e 

Thus, {6,x,y) are local coordinates on SE{2). 

A basis of the Lie algebra se(2) = of SE{2) is given by 






ei = I U U U I , 62 

and we have that 

[e, ei] = 62, [e, 62] = -61, [ei, 62] = 0. 
An element ^ e se(2) is of the form 

^ = W 6 + Di 61 + U2 62 

and the exponential map exp : 5e(2) ^ ^ SE{2) of SE{2) is given by 
, . , sin Lo cos ll> — 1 cos w — 1 sin lu 

exp(w,Dl,U2) = (w,Ul \-V2{ ),-^^l( ) +V2 ), it uj ^0, 

and 

exp(0, wi,W2) = (0,i;i,W2). 

Note that the restriction of this map to the open subset U =] — 7r,7r[xIR^ C 9 
se(2) is a diS'eomorphism onto the open subset exp(J7) of SE{2). 

A discretization of the Chaplygin sleigh may be constructed as follows: 

- The discrete Lagrangian Ld : SE{2) — > M is given by 

Ld{n) = ^Tr (nm^) - Tr (m), 

where J is the matrix: 

(J/2) + ma^ mab 
mob (J/2) + mfo^ 

ma mb 

(see [H]). 

- The vector subspace Dc of se(2) is 

T>c = span {e, ei} = { (w, wi, U2) S se(2) | ^2 = } . 

- The constraint submanifold Mc of SE{2) is 
Mc = exp([/nDc). (4.6) 

Thus, we have that 

Mc = {(6',x,?/) G 5£;(2) I -TT < 6* < TT, e 7^ 0,(1 -cos6')x- ysine* = 0} 
U{(0,a;,0) G S'£;(2) | a; G M}. 
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Figure 1 . Submanifold Mc 



From (4.6 1 it follows that / £ Mc and T/Mc = Dc- In fact, one may prove 
that 

X 

T(o,x,o)^c^ span {— + o ^ ' a~ 

oO \{o,x,o) 2 5y|(o,:.,o) dx\{o,x,o) 



for X eR. 

Now, the discrete Euler-Poincare-Suslov equations are: 

^{9i,xi,yi){Ld) - ~e{02,X2,V2){Ld) = 0, 
e{{9i,Xi,yi){Ld) - el{d2,X2,y2){Ld) = 0, 

and the condition {9k,Xk,yk) € Mc, with k E {1,2}. We rewrite these equations as 
the following system of difference equations: 

—am cos 9i — bm sin 9i + am \ f mx2 + am cos 02 
+mxi cos 01 + mj/i sin j \ —bmsm92~am 

am 2/1 cos 01 — am xi sin 01 — 6r7i xi cos 01 \ / amy2 — bmx2 
—bmyism9i + {a^m + b^m+j)sm9i J y +{a^m + b^m + J) s'mt 

together with the condition 

iOk,Xk,yk) e Mc, k(E {1,2}. 

On the other hand, one may prove that the discrete nonholonomic Lagrangian 
system {Ld,M.c,T>) is reversible. 

Finally, consider a point {0,x,0) E Mc and an element 77 = {uj,vi,V2) G se(2) 
such that 

V(0,x,0)er(o,,,o)Mc, V(0,a;,0)(l?(Ld)) = 0, ^(0, 0)(i^(Ld)) = 0. 

Then, if we assume that a^m + J + am^ 7^ it follows that rj — 0. 

Thus, the discrete nonholonomic Lagrangian system (L(i,Mc,2)c) is regular in a 
neighborhood of the identity /. 

4.4. Discrete nonholonomic Lagrangian systems on an action Lie group- 

oid. Let i7 be a Lie group with identity element e and ■ : M x H ^ AI, {x, h) S 
M X H 1-^ xh, a right action of H on M. Thus, we may consider the action Lie 
groupoid T = M X H over M with structural maps given by 

a{x,h) = x, (3{x,h) = xh, e{x) = {x,e), 

fh{{x, h), {xh, h')) — {x, hh'), i{x, h) — {xh, h~^). 
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Now, let t) — T^H be the Lie algebra of H and $ : f) ^ X(Af) the map given by 

= riM, for 77 G [), 

where rjM is the infinitesimal generator of the action • : M x H ^ M corresponding 
to rj. Then, $ is a Lie algebra morphism and the corresponding action Lie algebroid 
pri : M X I) M is just the Lie algebroid ofT = M x H. 

We have that Sec(pri) = {fj : M ^ t) \ fj is smooth } and that the Lie algebroid 
structure (|-, -J^jp^) on pri : M x H ^ M is defined by 

lfl,lLl^{x) = [f]{x),ii{x)]+{fi{x))M{x){jx)-{fi{x))M{x){fi), p^{fi){x) = {fj{x))M{x), 
for fi,flG Sec(pri) and x G M. Here, [•, •] denotes the Lie bracket of i). 

If {x,h) E T = M X H then the left-translation l(^x.h) ■ a^^ixh) — > a^^{x) and 
the right-translation r^^-jj) : l3^^{x) — s- I3~^{xh) are given 

l^^,^h){xKh')^{x,hh'), r^,^h){x{h')-\h') = {x{h')-\h'h). (4.8) 
Now, if ?7 g f) then 77 defines a constant section C,, : Af ^ f) of pri : M x i) ^ M 



and, using (2.4|, (2.5), (4.7l and (4.8l, we have that the left-invariant and the 
right-invariant vector fields and C\, respectively, on M x H are defined by 

C^{x, h) = {-ijM{x),lf{h)), C^{x, h) = (0,, (4.9) 

for {x, h) gT = M X H. 

Note that if {rji} is a basis of f) then {C\.} is a global basis of Sec(pri). 

On the other hand, we will denote by expp : Er = M x f) ^ F = Af x _ff the 
map given by 

expp(a;,?7) = (a;, expjj(ry)), for {x,r]) G Er = M x t), 

where expj:^ : i) H is the exponential map of the Lie group H. Note that if 
^(x,e) : R ^ r = Af X H is the integral curve of the left-invariant vector field C^j 



onT = M X H such that <i>(j._e)(0) = {x,e) then (see (4.9)) 

expr(a;,77) = <^(x,e){^)- 

Next, suppose that Ld : F = Af x iJ ^ M is a Lagrangian function, Dc is a 
constraint distribution such that {^"} is a local basis of sections of the annihilator 
T>^, and Mc C F is the discrete constraint submanifold. 

For every h G H (resp., x G M) we will denote by L/j (resp., L^) the real 
function on M (resp., on H) given by Lh{y) — Ld{y, h) (resp., Lx{h') = Ld{x, h')). 
A composable pair {{x, hk) , {xht, hk+i)) G F2 n (Mc x Mc) is a solution of the 
discrete nonholonomic Euler-Lagrange equations for the system (Lc(,Mc,'Dc) if 

Cjj{x,hk){Ld) - C.,j{xhk,hk+i){Ld) = XaX°'{xhk){r]) , for all 77 G f), 



or, in other terms (see (4.9)) 
for all 77 G i). 

4.4.1. The discrete Veselova system. As a concrete example of a nonholonomic 
system on a transformation Lie groupoid we consider a discretization of the Veselova 
system (see |11]). In the continuous theory [9], the configuration manifold is the 
transformation Lie algebroid pri : S"^ x so (3) with Lagrangian 

Lc(7, uj) ^ ^uj ■ luj - mglj ■ e, 
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where S'^ is the unit sphere in M.^, a; e M"^ ~ so(3) is the angular velocity, 7 is the 
direction opposite to the gravity and e is a unit vector in the direction from the 
fixed point to the center of mass, all them expressed in a frame fixed to the body. 
The constants m, g and I are respectively the mass of the body, the strength of 
the gravitational acceleration and the distance from the fixed point to the center 
of mass. The matrix I is the inertia tensor of the body. Moreover, the constraint 
subbundle Dc S"^ is given by 

7 e S'^ i-> Dc(7) = { w G ~ so(3) I 7 • w ^ } . 

Note that the section (j) : ^ S'^ x 5o(3)*, {x, y, z) i-> {{x, y, z),xe^ + ye^ + ze^), 
where {ei, 62, 63} is the canonical basis of and {e^, e^, e^} is the dual basis, is a 
global basis for 

If w £ so (3) and u is the skew-symmetric matrix of order 3 such that ujv — lu x v 
then the Lagrangian function L^, may be expressed as follows 

Lc{j, ^) = ^ Tr(a} JS^) -mglj-e, 

where I = ^ Tr(/)/3x3 — /. Here, /3X3 is the identity matrix. Thus, we may define 
a discrete Lagrangian Ld : T = S"^ x 5*0(3) R for the system by (see [57]) 

Ld(7,fi) = -~^Tr{m) - hmgl-^-e. 

On the other hand, we consider the open subset of 5*0(3) 

V={n€SO{?>) |Trf]^±l} 
and the real function ijj : S'^ x V ^ M. given by 

V'(7,^^) = 7•(^^'^^)■ 

One may check that the critical points of ijj are 

O^, = {(7,^^) e 52 X 1/ I 177-7 = 0}. 
Thus, the subset of T = 5^ x 50(3) defined by 

= { (7,17) G (52 X V^) - O^ I 7 • (IWF) = 0} , 

is a submanifold of F of codimension one. Mc is the discrete constraint submanifold. 

We have that the map expp : x so(3) x 50(3) is a diffeomorphism from 

an open subset of Dc, which contains the zero section, to an open subset of Mc, 
which contains the subset of T given by 

e(52) = {(7,e)G 5^x50(3)}. 

So, it follows that 

(Dc)(7) =T(^,e)Mcni;r(7), for 7 G 5^. 

Following the computations of [37] we get the nonholonomic discrete Euler-Lagrange 
equations, for {{-fk,^k), {-fk^k,^k+i)) G r2 

Mk+i-9^Mk^k + mglh'^{-ij^r^ e) = Xjk+i, 

7fe(rifc - = 0, 7fc+i(rifc+i - rife+i) = o, 

where M — ill — lil^ . Therefore, in terms of the axial vector H in M.^ defined by 
n = M, we can write the equations in the form 

Ilfe+i = ^^fc^fc - mglK^^k+i X e + A7fc+i, 

= 0, 7fc+i(f^fc+T^I^+i) = 0. 
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Note that, using the expression of an arbitrary element of 5*0(3) in terms of the 
Euler angles (see Chapter 15 of [31]), we deduce that the discrete constraint sub- 
manifold Mc is reversible, that is, i(Mc) = Mc- However, the discrete nonholonomic 
Lagrangian system (Lj_,JAc,Di.) is not reversible. In fact, it is easy to prove that 
Ldoi ^ Ld- 

On the other hand, if 7 G 5^ and ^, r] eR^ = so(3) then it follows that 
Ce(7,/3)(C,(L,))--e-/r;. 

Consequently, the nonholonomic system {Ld,M.c,'Sc) is regular in a neighborhood 
(in Mc) of the submanifold e{S^). 

4.5. Discrete nonholonomic Lagrangian systems on an Atiyah Lie group- 

oid. Let p : Q M ~ Q/G be a principal G-bundle and choose a local trivial- 
ization G x U, where U is an open subset of M. Then, one may identify the open 
subset ip-^{U) X p-\U))/G ~ ((G x [/) x (G x U))/G of the Atiyah groupoid 
{Q X Q)/G with the product manifold {U xU) x G. Indeed, it is easy to prove that 
the map 

((G X J7) X (G X U))/G -^{UxU)xG, 

[{{g,x),{g\y))]^{{x,y),g-'9')), 

is bijective. Thus, the restriction to ((G x J7) x (G x U))/G of the Lie groupoid 
structure on {Q x Q)/G induces a Lie groupoid structure in ([/ x [/) x G with 
source, target and identity section given by 

a : {U X U) X G ^ U; ((a;, y), g) — > x, 
f3:{UxU)xG^U; {{x,y),'g) ^ y, 
e:U^{UxU)xG; x^{{x,x),t), 

and with multiplication m : {{U x U) x G)2 {U x U) x G and inversion i : 
{U xU) X G ^ (U X U) X G defined by 

m{{{x,y),9)A{y,z),h)) = {{x,z),gh), 
iiix,y),g) = {{y,x),g-^). 

The Lie algebroid A{{U xU)xG) may be identified with the vector bundle TU x g — s- 
U . Thus, the fibre over the point a; G C/ is the vector space T^U x g. Therefore, a 
section of A{{U x [/) x G) is a pair (A, ^), where A is a vector field on U and ^ is a 
map from U on g. The space Sec{A{{U x U) x G)) is generated by sections of the 
form {X, 0) and (0, Gj), with X e X.{U), ^ € and G^ : [/ — > g being the constant 
map C^{x) — ^, for all x £ U (see [27] for more details). 

Now, suppose that : (f^ x C^) x G ^ M is a Lagrangian function, Dj. a vector 
subbundle of TU x g and Mc a constraint submanifold on {U xU) xG. Take a basis 
of sections {Y"} of the annihilator D°. Then, the discrete nonholonomic equations 
are 

< > 

{Xa,fia){{x,y),gk){Ld) - {Xa,f}o){{y, z),9k+i){Ld) = 0, 
with {Xa,f}a) ■ U TU X a basis of the space Sec(TD^) and {{{x,y),gk),{{y, z), 
9k+i)) G (Mc X Mc) n {{U X U) X G)2- The above equations may be also written as 



{X, 0){ {x, y),gk){Ld) - { X, 0)(( y, z),gk+i){Ld) = Xc.Y''{y){X{y)), 
(0, Ge)((x, y),gk){Ld) ~ (0, C^){{y, z),gk+i){Ld) - KY''{y){C^{y)), 

with X e X{U), e e and (((x, y),gk), {{y, z), fffe+i)) G (Mc x Mc) n {{U xU)x G^. 
An equivalent expression of these equations is 

D2Ld{{x,y),gk) + DiLd{{y, z), gt+i) = AaAi"(y), ,^ 
Pk+i{y,z) = Ad*g^Pk{x,y) - Xafi°'{y), 
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where pk{x,y) = d{r*^L(^x.y, ))(^) for {^iV) ^ U x U and we write Y" = (^",7)"), 
/Lt" being a 1-forni on U and fy" : [/ g* a smooth map. 

4.5.1. A discretization of the equations of motion of a rolling ball without sliding on 
a rotating table with constant angular velocity. A (homogeneous) sphere of radius 
r > 0, mass m and inertia about any axis / rolls without sliding on a horizontal 
table which rotates with constant angular velocity about a vertical axis through 
one of its points. Apart from the constant gravitational force, no other external 
forces are assumed to act on the sphere (see [H]). 

The configuration space for the continuous system is Q = x 50(3) and we shall 
use the notation (x, y; R) to represent a typical point in Q. Then, the nonholonomic 
constraints are 

x+^-Tt{RR^E2) = -ny, 

y-'^TriRR^Ei) = ^x, 

where {Ei, E2, E^} is the standard basis of 5o(3). 

The matrix RR^ is skew symmetric, therefore we may write 

/ -W3 W2 \ 
RRT ^ \ W3 -Wi 

\ -W2 Wi / 

where {wi,W2,W3) represents the angular velocity vector of the sphere measured 
with respect to the inertial frame. Then, we may rewrite the constraints in the 
usual form: 

X — rw2 = —^y, 
y + rwi — fix. 

The Lagrangian for the rolling ball is: 

L,ix,y;R,x,y;R) = ^m{x^ + f) + ^ITriRR^ (RR^f) 

= ^m{x^ + f) + ^I{u;j+u;l+ujl). 

Moreover, it is clear that Q = x S0{3) is the total space of a trivial princi- 
pal S'0(3)-bundle over and the bundle projection (/> : Q — > M = is just the 
canonical projection on the first factor. Therefore, we may consider the correspond- 
ing Atiyah algebroid E' ^ TQ/SOii) over M = identify the tangent 
bundle to 50(3) with so(3) x 50(3) by using right translation. 

Under this identification between r(50(3)) and so(3) x 50(3) the tangent action 
of 50(3) on T(50(3)) =^ 5o(3) x 50(3) is the trivial action 

(50(3) X 50(3)) X 50(3) ^ so(3) x 50(3), {{uj, R),S) k-^ {lu, RS). (4.12) 

Thus, the Atiyah algebroid TQ/ SO{S) is isomorphic to the product manifold 
TM^ X so (3) and the vector bundle projection is Tr2 opri , where pri : TM? x so (3) — ^ 
TM? and tr2 : TM? are the canonical projections. 

A section of E' = TQ/SO{3) ^ TM^ ^ so(3) ^ MMs a pair {X,u), where X is 
a vector field on and m : ^ so (3) is a smooth map. Therefore, a global basis 
of sections of TM^ x so (3) 

4 = (i),Bi), »; = (i),Ej), 4 = (o,£3). 
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The anchor map p' : E' ^ TQ/SO{S) ^ TM^ x so(3) TM^ is the projection 
over the first factor and if |-, •]' is the Lie bracket on the space Sec{E' = TQ/ SO{3)) 
then the only non-zero fundamental Lie brackets are 



4,41' 



'5, 



14,41' 



'3, 



14,41' 



Moreover, the Lagrangian function Lc = T and the constraint functions are 
S'0(3)-invariant. Consequently, Lc induces a Lagrangian function on E' = 
TQ/SO{i) 



where {x,y,x,y) are the standard coordinates on TM^ and uj G so(3). The con- 
straint functions defined on E' — TQ/ SO{3) are: 

i+§Tr(wS2) = -ny, 
2/ — I Tr(a;i?i) — fix. 



(4.13) 



We have a nonholonomic system on the Atiyah algebroid E' ~ TQ / SO{3) = TM^ x 
so(3). This kind of systems was recently analyzed by J. Cortes et al 9 (in particular, 
this example was carefully studied). 

Eqs. (4.131 define an affine subbundle of the vector bundle E' ^ TM? x so(3) — > 
M? which is modelled over the vector subbundle D'^ generated by the sections 

2)^ = {4,r.'i + 4,r4-4}. 

Our objective is to discretize this example directly on the Atiyah algebroid. The 



Atiyah groupoid is now identified to E"^ x x SO{i) 
the discrete Lagrangian by 



L'd{xo,y(i,xi,yi;Wi) = L'c{xo,ya 



xi -xo yi- yo 



We may construct 



(logWi)/h) 



h ' h 

where log : 50(3) — > so (3) is the (local)-inverse of the exponential map exp : 
so(3) — > 50(3). For simplicity instead of this procedure we use the following 
approximation: 

log Wi/h « 

where /3X3 is the identity matrix. 
Then 



L[i{xo,yQ,xi,yi;Wi) = L'c{xo,yo 



xi - Xo yi- yo Wi - /; 



3x3 



xi - Xo 



h ' 

yi - yo 



■ Tr(/: 



3x3 



Eliminating constants, we may consider as discrete Lagrangian 



Xl - Xo 

h 



yi - yo 
h 



2h^ 



Tr(VFi) 



The discrete constraint submanifold MJ, of M-^ x x 50(3) is determined 
by the constraints: 



xi - Xo 
h 

yi - yo 



2h 

r 



Tr{WiE2 



-n 



yi + yo 



2 

Xi + Xo 
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We have that the system {L'^,l\i'^,D'^) is not reversible. Note that the Lagrangian 
function is reversible. However, the constraint submanifold is not reversible. 

The discrete nonholonomic Euler-Lagrange equations for the system (L'^,,M'^, 
D^) are: 

< — — > 

< > 

{rs[ + S4)(a;o,yo,a;i,yi; VFi)(L^) - {rs[ + s'4^){xi,yi, X2,y2;W2){L'j) = 

< > 

{rs2~ s'^){xQ,yo,Xi,yi;Wi){L'^) - {rs2~ s'^){xi,yi,X2,y2;W2){L'^) = 
with the constraints defining JAc- 

^ ^ < > < 

On the other hand, the vector fields s g, s 5, rs'i + S4, rs'i + S4, rsj — S3 and 

rs'2 - s'^ on (M^ X X SO{3) are given by 

V^(0,0),^3), - ( (0,0), ^ 3), 

gT4 = ((0,rJ|),^2), r s[ + si = {{~r^,0),E2), 

rs'2 - 4 = mr^),-E,), rs'2 - 4 = m-r-§^),~E,), 

where Ei (respectively, Ei) is the left- invariant (respectively, right-invariant) vec- 
tor field on 5*0(3) induced by Ei £ so(3), for i € {1,2,3}. Thus, we deduce the 
following system of equations: 



Tr ((14^1-1^2)^3) - 0, 
+ ■ ' Tr{{W,-W2)E2) = 0, 

Tr ((14^1-1^2)^1) = 0, 



; 2/i2 

y2 - 2yi + yo \ I_ 

2h? 









X2 - 


x,[ 


r 


h 




2h 


2/2 - 




r 


h 




2h 



TT{W2E2) + n^^^^ = 0, 



2:2 -t- Xi 







where (xq, xi, yo, yi, Wi) are known. Simplifying we obtain the following system of 
equations: 

X2 - 2xi + xq ^ in y2 - yo 

/ + rnr^ 2h 

y2 - 2yi + yo _ in X2 - xo 

I + mr'^ 2h 

Tr ((11^1-14^2)^3) 

X2-X1 r ^2/2 + yi 

__ + _Tr(14^2i^2)+^^^^ 

y2-yi r ^x2 + xi 

^— --Tr(14^2i?l)-^^^^ 
Now, consider the open subset J7 of x ^ 5'0(3) 

U ^ {R^ xR^) X {W e 50(3) I W - Tr(14^)j 



= 


(4.14) 


= 


(4.15) 


= 


(4.16) 


= 0, 


(4.17) 


= 0. 


(4.18) 


x3 is regular}. 





Then, using Corollary 3.13 (iv), we deduce that the discrete nonholonomic La- 
grangian system {L'^,M'^,T>'^) is regular in the open subset U' of given by 

U' = f/nM' 
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If we denote by — {xk+i — Xk)/h and Vk — {Vk+i — yk)/h, fc e N then from 
Equations (4.141 and (4.15) we deduce that 



Uk+l 
Vk+l 

or in other terms 



= A 



Uk 
Vk 



1 



-Aah 



Uk 
Vk 



x{k + 2) = 
y{k + 2) = 



8x{k + 1) + {a^h^ - 4:)x{k) - Aah{y{k + 1) - y{k)) 
8y{k + 1) + (a^h^ - A)y{k) + Aa{x{k + 1) - x{k)) 



Since A e 50(2), the discrete nonholonomic model predicts that 



where a — , , 

the point of contact of the ball will sweep out a circle on the table in agreement 
with the continuous model. Figure [2] shows the excellent behaviour of the proposed 
numerical method 





Figure 2. Orbits for the discrete nonholonomic equations of mo- 
tion (left) and a standard numerical method (right) (initial condi- 
tions a;(0) = 0.99, y{0) = 1, x(l) = 1, y(l) ^ 0.99 and h = 0.01 
after 20000 steps). 



4.6. Discrete Chaplygin systems. Now, we present the theory for a particu- 
lar (but typical) example of discrete nonholonomic systems: discrete Chaplygin 
systems. This kind of systems was considered in the case of the pair groupoid in 

m- 

For any groupoid F =^ M, the map x : T ^ M x M, g i-^ {a{g), l3{g)) is a 
morphism over M from F to the pair groupoid M x M (usually called the anchor 
of F) . The induced morphism of Fie algebroids is precisely the anchor p : Er TM 
of Et (the Fie algebroid of F). 

Definition 4.2. A discrete Chaplygin system on the groupoid T is a discrete 
nonholonomic problem (i^j, Mc, Dc) such that 

- (Lrf, Mc, Dc) is a regular discrete nonholonomic Lagrangian system; 
~ XMc — X° *Mc • — > M X M is a diffeomorphism; 

- p o i-X)^ : T)c — > TM is an isomorphism of vector bundles. 

Denote by Ld : M x M — > M the discrete Fagrangian defined by Ld ^ Ld o 

iMa ° iXMa)^^- 

In the following, we want to express the dynamics on M x M, by finding relations 
between de dynamics defined by the nonholonomic system on F and M x M. 
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From our hypothesis, for any vector field Y G X{M) there exists a unique section 
X g Sec(TD^) such that p o i^j^ o X = Y. 

Now, using and ([2^, it follows that 

Tga(x{g)) = -r(a(5)) and Tg(3{X{g)) = Y{(3{g)) 

with some abuse of notation. In other words, 

T,x(X(i'0)(5)) = Y^^'°\a{g),(3{g)) and T,x(^^°^'H.9)) - Y(^^^\a{g), P{g)) 

for g e Mc, where 7x : T^r ^ Vf3 ®r l^a T*^^*^(M x M) = T{M x M) is the 
prolongation of the morphism x given by 

i7gx){Xg,Yg) = {{Tga){Xg),{Tg(3){Y,)), 
for g e r and {Xg, Yg) G T^F = K,/? ® Fga. 

Since xMc is a diffeomorphism, there exists a unique G TgMc (respectively, 
G TgMc) such that 

{TgXMjiX'g) = r(i^")(a(g),/3(.9)) = (-r(a(5)), 0^(,)) 

(respectively, (TgXMj(^^) = Y(°'^\a{g), P{g)) = (0^(3), r(/3(5)))) for all g G M,. 
Thus, 

G TgMc n £(5) -x'g = z'ge Vga n yg/3, 

G TgMc n V^ga, X{g)-X'g^Z'geVganVg(3, 

for all g G Mc- 

Now, if {g, /i) G Fa n (Mc x Mc) then 

X(5)(L,)-X(;i)(Ld) = X;(Ld) + ^;(id)-X,;(id)-^;,(id) 

= Y{a{g),P{g)){La) - Y {a{h), P{h)){Ld) 
+Z'g{La)-Z\,{Ld). 

Therefore, if we use the following notation 

{a{g),(3{g)) = (x, y), {a{h),P{h)) = (y, z) 



then 



X{g){Ld) - X{h){Ld) - Y{x,y){Ld)~Y{y,z){Ld) 

~F+{x,y) + F-{y,z). 

In conclusion, we have proved that (g, h) is a solution of the discrete nonholonomic 
Euler-Lagrange equations for the system (Ld,Mc,I'c) if and only if {{x,y), {y,z)) 
is a solution of the reduced equations 

Y{x, y){Ld) - Y{y, z){Ld) = F+{x, y) - Fy{y, z), Y G X(M). 

Note that the above equations are the standard forced discrete Euler-Lagrange 
equations (see [22]). 

4.6.1. The discrete two wheeled planar mobile robot. We now consider a discrete 
version of the two-wheeled planar mobile robot [8,9. The position and orientation 
of the robot is determined, with respect a fixed cartesian reference, by an element 
Q, = {9, X, y) G SE{2), that is, a matrix 




= ( sin 6* cos 9 y 
1 
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Moreover, the different positions of the two wheels are described by elements 
{(jfjip) G T^. Therefore, the configuration space is SE{2) x T^. The system is 
subjected to three nonholonomic constraints: one constraint induced by the condi- 
tion of no lateral sliding of the robot and the other two by the rolling conditions of 
both wheels. 

It is well known that this system is S'£'(2)-invariant and then the system may 
be described as a nonholonomic system on the Lie algebroid se(2) x TT^ (see 
[3]). In this case, the Lagrangian is 



^Tr(an + y'^^ + yV:'^ 



where 





ei + v'^ €2 = \ I and 



Here, m = toq + 2mi, where niQ is the mass of the robot without the two wheels, 
nil the mass of each wheel, J its the moment of inertia with respect to the vertical 
axis, J2 the axial moments of inertia of the wheels and / the distance between the 
center of mass of the robot and the intersection point of the horizontal symmetry 
axis of the robot and the horizontal line connecting the centers of the two wheels. 

The nonholonomic constraints are 

wi + f</>+f^/' = 0, 

= 0, (4.19) 

determining a submanifold M of se(2) x TT^, where R is the radius of the two 
wheels and 2c the lateral length of the robot. 

In order to discretize the above nonholonomic system, we consider the Atiyah 
groupoid r = SE{2) X (T^ x T^) z4 t\ The Lie algebroid of SE{2) x (T^ x T^) =t 
is TT^ X 5e(2) -> T2. Then: 

- The discrete Lagrangian Lj. : SE{2) x (T^ x T^) ^ M is given by: 

Ld{ilk,(l)k,'(pk,(f>k+i,^k+i) = 2W^r {{nk ~ I3x3)^{^k - hxaV) 

I Ji (Mkf , JiiMji)! 

where /3X3 is the identity matrix, A4>k = (/Jk+i ~ (t>k, ^i^k = V'fe+i - V'fc 
and 

(cos 9k — sin ( 
sin Ok 

We obtain that 
Ld = ^ {mxl + myl- 2lmoXk{l- cos 9k) 

+2J{1 - cos 9k) +2lmoyk sm9k) + I'h^^j^ + l^^^^J^- 

- The constraint vector subbundle of se(2) x TT^ is generated by the sections: 
R R d R R d 
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The continuous constraints of the two-wheeled planar robot are written in 
matrix form (see 4.191: 














We discretize the previous constraints using the exponential on SE(2) 



(see Section 4.3.2) and discretizing the velocities on the right hand side 

i(^A0fc-i|AVfc) \ 



cos(i|A0fc-i|A'0fc) sin(i|A0fc-i|A^;.) -c 
-sin(i|A0,~i|A^,) cos(i|A0,--i|A^,) 



— Ai/- 



(l-cos(i|A0fc-i|A^fc)) 



if A</>fe 7^ A^Jk and 




if A(/)fc = At/jk. 

Therefore, the constraint submanifold Mc is defined as 

2c 2c 

A<^fe + A^.fc . (R.. R 
-''A<h^.''''\2c^'^'~2c 



Xk 



Atpk 



Vk 



A(t>k + AV'fc 



1 



COS I — A(/>fe 
2c 



R 



Ai^k 



(4.20) 
(4.21) 

(4.22) 



A(j)k~Aijk V V2c 2c 

if A(t>k ^ AtPk and 6ife = 0,Xk = -RAcPk and yfe = if A^ifc = A^-fe. 

We have that the discrete nonholonomic system (L^jMcDc) is reversible. More- 
over, if er : T-^ ^ SE{2) x (T^ x T^) is the identity section of the Lie groupoid 
r = SE{2) X (T^ X T^) then it is clear that 

er(T2)=={/3>,3}x At2xt2 CM,. 

Here, At2xT2 is the diagonal in x T^. In addition, the system {Ld,M,c,T)c) is 
regular in a neighborhood U of the submanifold er(Tr^) = {^3x3} x A^axT^ in M,. 
Note that 

^ihx3,'Pi,i>u'Pui>i)^c n Eri(t>i,ilJi) = Dc(0i,'0i), 

for {(j)i,ipi) e T^, where Er ~ se(2) x TT^ is the Lie algebroid of the Lie groupoid 
r = SE{2) X (T2 X T^). 

On the other hand, it is easy to show that the system (L^, U, Dc) is a discrete 
Chaplygin system. 

The reduced Lagrangian on x is 

1 ^A^k + Aih^n t^A^ ^ 

R . . R 1 . (A0fc)2 1 {Aijk? 



-AiPk)) 



-J{1 - cos(— A0fc 
2c 



2c 



AV'fe))) 



/l2 



/l2 



if A^fe^AV'fe 



, ^ mi?2 (A(^fe)2 



2 ' /i2 

The discrete nonholonomic equations are 



(Ld) - ~st 



(il2:<t>2:i>2,<t>3,i>3) 



(Ld) 



^2|(n2</'2,-</'2,</'3,'</'3)(^'*'' 
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These equations in coordinates are: 



JR, 



'2Ji{(t>3 — 202 + (jJi) ~ lRmo(cos 62 + cos ^i) H (sin 62 — sin ^i) 

c 



Rcost 



-{Imoyi + cmxi) + 



R sin 01 



R 

H — {cmx2 + lmQ{y2 - 2c)) 



{ImoXi — cmyi) 



JR, 



(4.23) 



2Ji('03 ~ 2^2 + ^1) = Zi?mo(cos6'2 +cos0i) (sin 6*2 

c 



sm f/i ] 



R cos 6*1 



i? 



{Imoyi — cmxi 



R sin 9i 



{IniQXi + cmyi) 



H — {cmx2 — l'mo{y2 + 2c)) 



(4.24) 



Substituting constraints ( |4.20| , ( |4.21[ ) and ( |4.22| in Equations ( |4.23| ) and ( |4.24[ ) 
we obtain a set of equations of the type = /i(0i, ^2, '/>3, V"!; "02, "03) and = 
/i(0i, 021 03i 01; V'2j V'3) which are the reduced equations of the Chaplygin system. 



5. Conclusions and Future Work 



In this paper we have elucidated the geometrical framework for nonholonomic 
discrete Mechanics on Lie groupoids. We have proposed discrete nonholonomic 
equations that are general enough to produce practical integrators for continuous 
nonholonomic systems (reduced or not). The geometric properties related with 
these equations have been completely studied and the applicability of these devel- 
opments has been stated in several interesting examples. 

Of course, much work remains to be done to clarify the nature of discrete non- 
holonomic mechanics. Many of this future work was stated in |36j and, in particular, 
we emphasize: 

- a complete backward error analysis which explain the very good energy 
behavior showed in examples or the preservation of a discrete energy (see 

El); , 

- related with the previous question, the construction of a discrete exact 
model for a continuous nonholonomic system (see [T7 l [32 | [36]): 

- to study discrete nonholonomic systems which preserve a volume form on 
the constraint surface mimicking the continuous case (see, for instance, 
[ini US] for this last case) ; 

- to analyze the discrete hamiltonian framework and the construction of 
integrators depending on different discretizations; 

- and the construction of a discrete nonholonomic connection in the case of 
Atiyah groupoids (see pilET] ). 

Related with some of the previous questions, in the conclusions of the paper of R. 
McLachlan and M. Perlmutter p5F, the authors raise the question of the possibility 
of the definition of generalized constraint forces dependent on all the points qk~i, 
Qk and Qk+i (instead of just qk) for the case of the pair groupoid. We think that 
the discrete nonholonomic Euler-Lagrange equations can be generalized to consider 
this case of general constraint forces that, moreover, are closest to the continuous 
model (see [^[55]). 
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